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ABSTRACT 

We consider the details of the QED processes that create electron-positron pairs in magnetic fields approach- 
ing and exceeding 10 14 G. The formation of free and bound pairs is addressed, and the importance of positro- 
nium dissociation by thermal X-rays is noted. We calculate the collision cross section between an X-ray and a 
| gamma ray, and point out a resonance in the cross section when the gamma ray is close to the threshold for pair 

conversion. We also discuss how the pair creation rate in the open-field circuit and the outer magnetosphere can 
be strongly enhanced by instabilities near the light cylinder. When the current has a strong fluctuating compo- 
nent, a cascade develops. We examine the details of particle heating, and show that a high rate of pair creation 
can be sustained close to the star, but only if the spin period is shorter than several seconds. The dissipation 
ON ' rate in this turbulent state can easily accommodate the observed radio output of the transient radio-emitting 

magnetars, and even their infrared emission. Finally, we outline how a very high rate of pair creation on the 
open magnetic field lines can help to stabilize a static twist in the closed magnetosphere and to regulate the loss 
' pH ' ' of magnetic helicity by reconnection at the light cylinder. 

O |. Subject headings: plasmas — radiation mechanisms: non-thermal — stars: magnetic fields, neutron 
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<M . 1. INTRODUCTION 

cza , 

Even when not emitting bright bursts of X-rays, magnetars continue to show dramatic variability on much longer timescales. 
Much of this variability appears to be a consequence of a dynamic magnetosphere. Various manifestations include extremely 
bright and persistent hard X-ray emission, which can outstrip the loss of rotational energy by a thousand times (Kuiper et al. 
2004, 2006); extreme changes in spindown torque that grow and decay over a period of months and can, in some cases, be 
^ I • sustained for years (Kaspi et al. 2001; Woods et al. 2002, 2007; Camilo et al. 2007a, Camilo et al. 2008); and, recently, the 
' discovery of high-frequency radio pulsations from two magnetars (Camilo et al. 2006; 2007b). A general review of magnetars 
jy-j \ can be found in Woods & Thompson (2006). 

^vq ■ This long-term variability appears to be connected only indirectly to the sudden, disruptive events known as Soft Gamma 
Repeater (SGR) bursts. Increases in spindown torque in the SGRs and some Anomalous X-ray Pulsars (AXPs) typically follow 
the X-ray bursts by a period of months. Very large torque variations are possible even if the transient X-ray output of the star 
is quite modest, as in the radio magnetar XTE J1810-197 (Camilo et al. 2007a). Even more remarkably, given the presence of 
strong torque variability, is the strong clustering of the spin periods of most magnetars within the range of 5 - 12 seconds. The 
decay of an internal magnetic field, which ultimately powers the magnetospheric activity, should not be sensitive to the details of 
current flow and pair creation in the magnetosphere. 
IT^i ' Basic energetic considerations suggest that most of the dissipation responsible for the very bright, hard X-ray emission of 
, magnetars is probably concentrated in the closed magnetosphere (Thompson & Beloborodov 2005; Baring & Harding 2007). The 
■ inner magnetosphere of a magnetar can sustain much stronger electric currents than the open-field circuit. Thompson, Lyutikov, 
& Kulkarni (2002) pointed out that starquakes create magnetospheric twists that store an enormous amount of energy in the 
toroidal magnetic field, enough to supply the observed non-thermal output for a period of years. Beloborodov & Thompson 
(2007) showed that the voltage in the inner magnetosphere is regulated by a continual e ± discharge. The electric field that 
accelerates the particles in the magnetosphere is sustained by the self-induction of the circuit, and the loss of energy by radiation 
induces a gradual decay in the twist. 

It has also been suggested that relatively strong currents flow within the outer magnetosphere of a magnetar, being sustained 
by the outward transfer of magnetic helicity from an inner zone where it is injected. The redistribution of the twist into the 
outer magnetosphere causes the poloidal field lines to expand slightly, thereby increasing the open-field voltage (Thompson et al. 
2002). This effect depends on the extraction of a tiny fraction the helicity that is stored in the inner magnetosphere, and therefore 
does not require a large net input of toroidal field energy into the magnetosphere. 

Considerations of the stability of a twisted magnetosphere, and observations of strong variations in the brightness of the radio 
magnetars on short timescales (Camilo et al. 2007a, 2008), suggest that the twist in the outer magnetosphere is strongly dynamic. 
A first investigation of pair creation in a dynamic magnetosphere was given by Lyutikov & Thompson (2005), with application 
to the radio eclipses of the double pulsar system J0737-3039. They noted that strong heating of the particles in the outer 
magnetosphere would result from a cascade to high spatial frequencies. 

In this paper we consider the formation of electron-positron discharges in the case where the outer magnetosphere is strongly 
dynamic. As a first step, we return to first principles and reconsider the basic mechanisms of pair creation in some detail. Previous 
theoretical work on e^ creation and radio death line for magnetars was inconclusive. Baring & Harding (1998, 2001) suggested 
that high-energy photons cannot create pairs because they split in the ultrastrong field. This suppression of pair creation would, 
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however, take place only if both polarization modes are able to split, a selection rule which is not supported by microphysical 
calculations of this process (Adler 1971; Usov 2002). In particular, resonantly scattered photons convert very quickly to e ± when 
B>4B Q . 

Magnetars are bright sources of thermal (keV) X-rays, and gamma rays are created at a high rate close to the star through 
resonant scattering of the X-rays. The scattered gamma rays propagating through the curved magnetic field can be expected 
to convert adiabatically to bound electron-positron pairs (Herold, Ruder, & Wunner 1985; Shabad & Usov 1986). We review 
this process and identify some peculiar features in super-QED magnetic fields. The dissociation of bound pairs depends on the 
presence of a source of ionizing photons (Usov & Melrose 1996; Potekhin & Pavlov 1997). In addition to photoionization by 
infrared photons, we identify the resonant scattering of thermal X-rays by one of the bound charges as an effective dissociation 
channel. 

We also reconsider the creation of pairs by collisions between gamma rays and thermal X-rays. The cross section for this pro- 
cess has been calculated previously using wave functions for the magnetized electrons and positrons that are not good eigenstates 
of the spin (Kozlenkov & Mitrofanov 1986). We recalculate the cross section, finding a slightly different result when the two 
colliding photons are both below the threshold for single-photon pair creation. We also identify a near resonance in the collision 
cross section when one of the photons is close to the pair creation threshold. 

The period clustering of magnetars might, at first sight, appear to mark out the location of a radio death line in the P — P 
plane. The combined results of this paper, and a companion paper (Thompson 2008), suggest the existence of at least two distinct 
regimes of pair creation: a more active regime in which the outer magnetosphere is a dynamic reservoir of magnetic helicity and 
the pair density is extremely high; and a more quiescent regime in which pair creation occurs at more modest rate, and the radio 
emission is concentrated closer to the star. In this paper, the mechanism of particle heating in the dynamic state is considered 
self-consistently, along with the energy and multiplicity of the heated pairs. Pair creation occurs mainly via resonant scattering 
of thermal X-rays close to the surface of the star. A critical spin period arises from the requirement that the heated charges can 
emit enough pair-creating gamma rays to compensate their loss from the circuit. 

Many features of the radio magnetars - their strong flux variability, hard radio spectra, and broad pulses - connect them with 
the more dynamic circuit behavior described here. The voltage that can be achieved on the more quiescent circuit is examined 
in detail in Thompson (2008), where it is shown that the power dissipated can barely explain the observed output of the radio 
magnetars, but is quite sufficient for more standard pulsar emission. The spin clustering of magnetars then signifies a limiting 
spin period for the more dynamic and intense regime of pair creation. Pulsed radio emission in the more quiescent regime should 
be possible at significantly longer spin periods (e.g. Medin & Lai 2007). 

The plan of the paper is as follows. We consider the conversion of resonantly scattered photons to pairs in and collisions 
between gamma rays and ther mal X -rays in $3] The formation of pairs on a magnetic flux rope with a dynamic twist is treated in 
$U An explanation is given in 35.21 for which current driven instabilities in the outer magnetosphere could be regulated by a high 
rate of pair creation on the open magnetic field lines. Possible sources of short-timescale variations in the open-field voltage are 
examined in $5] It is shown in $6]that most of the accumulated spindown of a magnetar can occur when its magnetosphere is in 
the active state. The paper closes with a summary of our results in $7] The Appendix provides details of the calculation of the 
matrix elements for 7 + 7 — > e + + e~ in an ultrastrong magnetic field. 

2. CREATION OF PAIRS VIA RESONANT SCATTERING OF THERMAL X-RAYS 

The voltage on the open magnetic field lines of a neutron star (and in the closed magnetosphere if it is non-potential) is regulated 
by pair creation. We focus in this section and $3] on the creation of pairs by energetic gamma rays. The background magnetic 
field will generally be assumed to be comparable to, or stronger than, the QED field 

Bq = V H^L_ — 4 4 x G. (1) 
eTi 

Magnetars are bright thermal X-ray sources, which means that scattering of the thermal X-rays is the dominant source of pair- 
creating gamma rays. This process has been discussed extensively in the context of radio pulsars, by Kardashev et al. (1984), 
Daugherty & Harding (1989), Chang (1995), Sturner 1995), and Hibschman & Arons (2001a,b). Macroscopic electric fields are 
screened rapidly enough that the electrons and positrons in the inner magnetosphere do not reach the energies that are required 
for gamma ray emission by curvature radiation. The electrostatic gap model investigated by Hibschman & Arons (2001a,b) and 
by Thompson (2008) has this property, as does the more dynamic circuit model described in this paper. 

We assume that the target X-rays are approximately described by a Planck spectrum with temperature 7bb, that is, by a spectrum 
with an exponential cutoff at Hujx 3> 3^^. This assumption is consistent with measurements of the X-ray spectrum of XTE 
J18 10— 197 (Halpern & Gotthelf 2005; Gotthelf & Halpern 2005). The spectra of other, brighter, magnetars have a significant tail 
at 1-10 keV, well above the typical k^Thb = 0.4 keV. However, this tail is likely generated at larger radii (Thompson et al. 2002; 
Lyutikov & Gavriil 2006; Fernandez & Thompson 2007) and diluted near the star compared with the thermal radiation from the 
star surface. 

Denote the angle of photon propagation with respect to B by 8^ and let /1 = cos^b- For definiteness, let us assume that the 
magnetic field is predominantly radial near the star. The stellar radiation at a radius r occupies the range of angles /i m ; n < /i < 1 
where 

/ R 2 \V2 

Mmi„W=fl--^J . (2) 

and 7?ns is the stellar radius. We neglect here the effects of gravitational light bending. The energy Viuj of a target photon is 
Doppler-boosted in the rest frame of the electron to 

hw' =7(l-/3fj,)hw. (3) 
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In this frame, the target photon is strongly aberrated and moves nearly parallel to B. 

The scattering cross section is strongly enhanced when the electron is promoted to its first Landau state. The resonant fre- 
quency 1 in the electron rest frame is (e.g. Herold 1979) 

, eB 

J = . (4) 

m e c 

This translates into a condition on the energy of target photons in the lab frame, which we write as 

_ _ huj b/b q 



kuTbb 6 b b 7(1 -/?/•*) ' 

where 



(5) 



Hbb = j (o) 

m e c z 

is the dimensional temperature of the stellar radiation. Setting f3 = 1 is a good approximation if the resonant scattering occurs 
close to the star, where HeB/m e c ^> ^B^bb- 

The target photons move nearly parallel to B in the rest frame of the scattering particle. For the cross section calculation, 
one can adopt the description of resonant scattering as absorption followed by de-excitation. Then the cross section is given by 
(Daugherty & Ventura 1978), 

/ r I I \ s 



7(1 - M ) 



(7) 



where uo B = eB/m e c and 



°L = ~-^- (8) 



27r 2 e 



The magnitude of the scattering cross section integrated through the resonance is 

~ ~ |^ (£) ^»a T , (9) 

where a/ ~ 1/137 is the fine structure constant. Note that <r res is much larger than the Thomson cross section oj even in a ~ 10 15 
G magnetic field. By comparison, the cross section below the resonance is strongly suppressed, 



a j 

"32' 



(10) 
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The two expressions here apply to both X-ray polarization eigenmodes. These modes are linearly polarized close to the star, 
where vacuum polarization effects dominate the dielectric properties of the medium. 

The rate of resonant scattering is high over a broad range of particle kinetic energies, due to the broad band nature of the black 
body radiation field. An electron moving radially away from the star with a Lorentz factor 7> ^ 1 scatters the thermal X-rays at 
the rate ^ 

r r >= J duj J 27rc/^^-(l-/i)cr res (w,M,7>), (H) 



where er res is the scattering cross section, and I u is the spectral intensity of thermal X-rays, 

Mb) 3 

87T 3 C 2 /i 



~ ^ 3 ^_ hu> 

/(D) = - s — , u> = — — . (13) 

We count only one polarization mode in I u as is appropriate to thermal emission from a strongly magnetized neutron star (e.g. 
Silantev & Iakovlev 1980). 

In the case of an outgoing electron near the neutron star surface, the double integral in equation dTTb then yields (S turner 1995) 



r res_ a / Q bb n .-,2-2 



2(B/Bq) 

where 



ln(l-e- Wc ") 



(14) 



w>(7>) = 7— t^ = — ^— 7z r, (15) 

KB Ibb 7> C'bb( 1 - Mmin) 

1 The resonant frequency is given by the non-relativistic cyclotron formula even in super-QED magnetic fields, because the absorbing charge experiences a 
recoil along B. Note also that a photon incident along the magnetic field can be resonantly scattered only at a single frequency, corresponding to the transition 
between the lowest and the first Landau states. Resonant interactions involving higher Landau states drop out for zero inclination angle. 
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and 7i> is the Lorentz factor of the scattering charge, aj = e 2 /He is the fine structure constant, and A7 = fi/m e c. The electron 
resonantly scatters soft X-rays from the Rayleigh- Jeans tail of the spectrum when its Lorentz factor is higher than 



7^ 



B/B Q 



(l _ Mmin)@bb 

The corresponding expression for an ingoing charge is 



= 2.3 x 10 4 



5, 



^B^bb 

0.5 keV 



where 



w< = 



__B/Bo_ 

fcfi7bb 29bb7re: 



= 1.2 X 10 4 



B] 



fa?bb N ' 
0.5 keV / 



(16) 



(17) 



(18) 



and we have made the approximation r^>R 



NS- 



2.1. Direct Conversion of Gamma-Rays to Pairs 

Once created in an ultra-strong magnetic field, a gamma ray will rapidly convert to e^ pairs. The details of this process 
depend sensitively on the strength of the magnetic field. The interval between gamma ray emission and pair creation is shortened 
substantially if the gamma ray is created above the threshold energy for pair creation. Photons this energetic are created by 
resonant scattering if the magnetic field is stronger than a threshold value, and are always created at some finite rate by non- 
resonant scattering. We call this process "direct" pair creation to distinguish it from the adiabatic conversion process. 

Resonant scattering may be viewed as an excitation of particle to the first Landau level followed immediately by de-excitation. 
The energy of the scattered photon depends on its pitch angle with respect to the magnetic field, . It is convenient to view this 
process in the frame where the intermediate (excited) electron does not move along B. In this frame, the energy of scattered 
photon is given by (Beloborodov & Thompson 2007), 



£ 7 = ■ 



sin 



1- 



cos + ■ 



2 4 

mic 



sin 



1/2" 



where Eg is the energy of the first Landau level, 



E B =[ 1+2 



B< 



1/2 



m e c 2 . 



The scattered photon will directly convert to an e pair if its energy exceeds the threshold 



E y > E t h T = 



2m e c 



(19) 



(20) 



(21) 



Photons scattered with = tt/2 have the minimum Etbr = 2m e c 2 and the maximum E 1 = Eb~ m e c 2 (because the de-exciting electron 
experiences no recoil). The condition £ 7 > is satisfied in a range of angles if EB-m e c 2 > 2m e c 2 , which is equivalent to 
B > 4Bq. If the surface magnetic field is stronger than this, a large fraction of resonantly scattered photons will directly convert 
to e ± near the surface. There is, of course, no constraint on the magnetic field for direct pair creation via non-resonant scattering. 

2.2. Adiabatic Conversion of a Single Gamma Ray to an Electron-Positron Pair 

A gamma-ray that is formed with an energy ftu> 3> 2m e c 2 and pitch angle O^b <S 2m e c 2 /Hu> is initially below threshold for pair 
creation, but can propagate through the curved magnetic field to a position where the condition (l2Tb is satisfied (Sturrock 1971). 
This process is familiar from models of ordinary pulsars with fields B^s = 10 12 - 10 13 G. Pair creation by a single photon can 
occur just above the threshold energy ( f2Lb only if the magnetic field exceeds a minimum strength, 



B 



4 B, 



Q 



: 0.05 ft 



Q 



:2x 10 12 G. 



(22) 



(e.g. Erber 1966; Berestetskii et al. 
state is 



3 ln[0.2a em 9 kB (rm e c/H)] 
1982). The corresponding maximum radius for the creation of a pair in the lowest Landau 



8fi 



1/3 

NS.15- 



(23) 



Pair creation by energetic photons, E 7 ^> 2m e c 2 , can occur above threshold in weaker magnetic fields, which leads to the emission 
of multiple synchrotron photons. 

In somewhat stronger magnetic fields, the gamma-ray is energetically favored to convert to bound positronium. We highlight 
here some ambiguities in the existing treatments of this process (Wunner et al 1985; Shabad & Usov 1986) when B > Bq. The 
associated refractive effects have been studied by Shabad & Usov (1984), who argue that the group velocity of the gamma ray 
will bend and asymptote to the local magnetic field. The perpendicular energy Tllo sin 9kB would, in the process, become very 
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close to 2m e c 2 . Qualitatively, a gamma ray that enters the dispersive regime at hk± ~ 2m e c can be viewed as a mixture of a pure 
electromagnetic mode (a photon) and a second mode (an electron-positron pair) that does not propagate across the magnetic field. 

The dispersion relation of the gamma ray receives a divergent correction from vacuum polarization near the threshold for pair 
creation. Neglecting the interaction between the electron and positron, one finds (Shabad & Usov 1984) 



ui 2 = c 2 k 2 -2a f I — I exp 



B 



-(kj_c/uj m ) 



21 



2(B/B Q ) 



(4w2- a; 2 + c 2 fe 2 ) l/2' 



(24) 



where uj„, = m e c 2 /H = 2nc/\. This dispersion relation applies to the two polarization modes which are luminal at low and high 
frequencies, and correspond to a propagating O-mode photon. The lower branch bends over from u> ~ ck at small momentum, to 
oj ~ m e c 2 /h when k± > m e c/h. 
The electrostatic force between the electron and positron will reduce the energy threshold for pair creation, 

2m e c 2 -E bind (k ± ,B) 

£W- -> — a • (25) 

sin 6^ 

The binding energy Ebmd of a positronium atom deviates significantly from the zero-field value i x 13.6 eV when B > fi a tomic = 
al m B Q = 2.3 x 10 9 G, 
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£bind — ■KZO* L z m m e c L ln z 

96 



B/B a 



l+(B/B Q )- l (hk ± /m e c) 2 



(26) 



(Wunner & Herold 1979; Shabad & Usov 1986; Lieb et al. 1992). The effect of this interaction on the gamma ray dispersion 
curve have been analyzed by Wunner et al. (1985) and Shabad & Usov (1985) by calculating the vacuum polarization using the 
positive energy wave functions of the bound electron and positron in the magnetic field. 

There are, however, some ambiguities in this procedure when B > Bq. First, the dispersion curve (l24l deviates strongly from 
that of a pure photon at the wavenumber k± = 2m e c/H. Writing hui = 2m e c 2 - AE on the lower branch, one finds 



AE If B 
2af B~Q, 



2/3 


4 


J exp 


. 3(B/B Q ) 



(27) 



One notices that AE has a strong exponential dependence on B, whereas £bind varies only logarithmically. In an ordinary pulsar- 
strength magnetic field, e.g. B = 4x 10 12 G, equation d27l i gives AE = 6 x 10~ 4 eV, much smaller than the positronium binding 
energy of ~ 100 eV. The electron-positron interaction is clearly important in this case, and the adiabatic conversion to a bound 
positronium atom is very plausible. In a super-QED magnetic field, the energy shift AE is much larger than the positronium 
binding energy as calculated from eq. 1261 : one finds AE = 15 keV when B = 4Bq. The energy levels of bound positronium, as 
calculated in the Bethe-Salpeter approximation, do not show the same distortion when hk± ~ 2m e c (Shabad & Usov 1986). 

Second, the propagation of a positronium atom in a curving magnetic field is qualitatively different from that of a photon. 
This process is allowed kinematically because the generalized momentum of each charged particle contains a term eA/c, where 
A = ^B x r is the background vector potential and r is the displacement from the center of mass. The conservation of momentum 
perpendicular to B gives 

1 2 

Ski-^^BxrJ, (28) 

a=l 

where a labels the electron and positron (<?i = — e, ei = +e). 

As the atom moves away from the star, the electron and positron are individually guided by the magnetic field and the perpen- 
dicular momentum of the atom decreases: equation (l28l implies k±^ ~ eBAr±^ oc B 1 / 2 . By contrast, the perpendicular wavevector 
of a photon generally increases as the photon propagates through an inhomogeneous magnetic field: the geometrical optics 
equations imply that k is essentially constant (even as the group velocity of the photon becomes aligned with the magnetic field). 

And, third, a direct calculation of the coulomb correction to the vacuum polarization gives the same dependence on B (Padden 
1994) only if B < B Q . 

In the absence of sufficiently high flux of target X-rays, the conversion of a gamma ray to bound positronium remains plausible, 
because AE decreases rapidly as k± rises above 2m e c/H. Once formed, a positronium atom can be dissociated by the resonant 
scattering of an ambient X-ray at a Landau resonance of the electron or positron ( jj2.3l ): or by the photoelectric absorption of an 
ambient IR photon ( 92.4) . A positronium atom in its ground state is not able to decay into two photons for a simple kinematic 
reason (Shabad & Usov 1982). The photon dispersion curve oj(k±) flattens near k± = 2m e c/h, which means that the conservation 
of momentum perpendicular to B is not consistent with 2 conservation of energy: two daughter photons with a total energy 
2m e c 2 - £bind must have a smaller momentum than was contained in the atom. 

2.3. Disintegration of a Positronium Atom through 
Absorption of an X-ray at the Electron/Positron Landau Resonance 

Given the ambiguities just mentioned in the dispersive properties of the hybrid photon-pair, we approach the problem of pair 
formation in two ways. The hybrid particle is treated either as a pure photon, or a pure positronium atom, and the relevant pair 



2 In other words, photon splitting is not facilitated by strong dispersion near the threshold for pair creation. 
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creation channels are analyzed. A gamma ray can collide with an ambient X-ray while it is still in the mixed state. We find that 
the cross section for such a two photon collision is resonantly enhanced when the gamma ray by itself is close to the threshold 
for pair conversion (50 ■ 

Although a positronium atom is electrically neutral, the electron and positron can scatter high-frequency photons whose wave- 
length is smaller than the electron-positron separation. The minimum separation is set by the transverse momentum of the atom: 
Arj_ ~ M±c/eB. A single electron at rest will resonantly absorb a photon of energy (B/Bq)m e c 2 with a cross section ||7). The 
corresponding photon wavelength A is 

Ar± 2 i y 2m e c} ^ 

We can therefore expect an X-ray photon to be absorbed by either the electron or positron with a cross section close to (|7). The 
photon imparts a recoil momentum to the absorbing particle, which is Ap ~ (B/BQ)m e c when B < Bq. The atom is dissociated 
if (Ap) 2 /2m e > Zibind, which requires 8>lx 10 12 G. Because bound pair creation occurs only for magnetic fields stronger than 
B > 2 x 10 12 G ~ 0.05 Bq, dissociation can be assumed to follow immediately following the resonant absorption of the X-ray. 

A positronium atom therefore has a very short mean free path for dissociation close to an X-ray bright magnetar. If the original 
gamma ray is created by resonant scattering of an X-ray by an electron of Lorentz factor 7, then the gamma ray has an energy 
— /recoil jm e c 2 , where the recoil factor is / reC oii = 5 when B = 4-Bq and / reco n ~ B/Bq for B <C Bq. The resultant positronium atom 
has a Lorentz factor r y± ~ (/recoil /2)7- The original seed photon for scattering, and the photon that dissociates the positronium 
atom, therefore tend to be drawn from different parts of the thermal X-ray peak. When B>Bq, the photons which act as seeds for 
pair creation come from above the black body peak, the dissociating photon comes from close to the peak, and the positronium 
is rapidly dissociated. On the other hand, when B < Bq, the dissociating photon is drawn from the Rayleigh- Jeans tail and, if 
£B7bb is below a critical value, it is possible for a newly created positronium atom to avoid dissociation by resonant scattering. In 
this circumstance the circuit voltage could be regulated by a different pair-creation channel - such as the direction conversion of 
non-resonantly scattered photons to pairs (Thompson 2008) - even though the creation rate of gamma rays by resonant scattering 
were formally higher. 

2.4. Photodisintegration of Parapositronium 

Many magnetars are intense sources of optical-IR photons, with a characteristic luminosity of 10 32 ergs s" 1 (Durant & van 
Kerkwijk 2006a,b, and references therein). These photons are in the right frequency range to dissociate bound positronium 
(parapositronium) when the relativistic motion of the bound pairs along B is taken into account. The binding energy is Eq ~ 
200 eV in a magnetic field B ~ 4Bq when the threshold for pair creation is Hk± = 2m e c/h (eq. l26l ). Pair-creating gamma 
rays typically have energies Tiuj ~ (0.05- \)^m e c 2 near the surface of the star, depending on the details of the circuit model 
(Thompson 2008). The resulting pair has a Lorentz factor 7 = hui/2m e c 2 ~ (500- 10 4 )Z?i5(A:B7bb/0.5 keV)" 1 parallel to B (see 
eq. |[T6l ). Photons of an energy 



fkj ~ - ~ (0.02 - 0.4)B7 eV (30) 

7(1-Mmin) V°- 5keV / 

will be most effective at photodissociating positronium. 

The cross section for photodissociation is suppressed compared with an unmagnetized positronium atom (see, e.g., Fig. 1 of 
Potekhin & Pavlov 1997). The photon is incident along B in the rest frame of the atom, and so both (linear) polarization modes 
have nearly identical cross sections. However, the motion of the charges in the atom is restricted in the plane perpendicular to B, 
and the cross-sectional area of the atom is reduced by a factor 



eB\ ( B 



- 



Tic j \B 



B V 00 .^-6^ B 



Here ag = 2h /m e e- is the Bohr radius of an unmagnetized positronium atom, and 

Ba t om 1 c=^a 2 m B Q = 6xl0 8 G (32) 

is the magnetic field that begins to significantly deform the atom away from its unmagnetized structure. 
The cross section for a photon incident parallel to B with an energy just above the threshold Eq can be estimated as 

7r 2 e 2 / B \~ l 
CTphot(w') ~ 7 7^77; /t , Tj » (33) 

{m e /2)(E a /K)C V^atomic/ 

where m e /2 is the reduced mass. (See Potekhin & Pavlov 1997 for more detailed expressions.) The rate of photodissociation in 
response to a spectral intensity I u of optical-IR photons is 

r p hot = J dui J 27rc/^(l-/i)^-CTphot(^'), (34) 
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as measured in the frame of the star. 

We now show that photodissociation is rapid if the optical-IR continuum of a magnetar is emitted close to its surface. In 
evaluating the above integral, we also make use of the observation that the luminosity density L u is roughly independent of 
frequency in at least some AXP sources. Then 

B ^ 



r p hot = (l-Mmin) -Tr „ ' ( 35 ) 

4r z m e c\E \ \B atomic J 

where L^/Anr 2 = 2^1^(1 - /i m i n ). It is convenient to express r p h t in terms of the luminosity ujL^ at Tioj = 0. 1 eV. Approximating 
1 ~~ Mmin — 5( r /^Ns)~ 2 ^ an d re-expressing r in terms of the local value of the magnetic field, gives 

1/3 

s" 1 . (36) 



7 (ojL u h 2 (By 



We see that photodissociation is rapid, r p h t > c/R^s, if the surface luminosity at 0.1 eV is higher than ~ 10 29 ergs s This 
critical luminosity is only ~ 10~ 2 - 10~ 3 of the typical K-band luminosities of AXPs; but it is also ~ 10 6 times larger than expected 
from the Rayleigh-Jeans tail of the X-ray black body. Photodissociation of positronium could therefore be slow if the optical-IR 
output were concentrated at a greater distance from the star. For example, a very dense pair gas forms on the open magnetic 
field lines in the circuit model described in $4] The pair density is high enough to seed the optical-IR emission of the transient 
magnetar XTE J1810— 197 by a form of plasma emission (Thompson 2008). In that case, the optical-IR flux would be beamed 
away from the star, and from the base of the open magnetic flux tube. 

3. COLLISIONS BETWEEN GAMMA RAYS AND THERMAL X-RAYS 

High-energy photons can create pairs by colliding with thermal X-rays. This pair creation channel is important when the X-ray 
flux is high (e.g. Zhang & Qiao 1998). Although most gamma rays will first convert to pairs off the magnetic field, a significant 
amount of pair creation by photon collisions may take place within a gap. If the surface temperature of the neutron star is higher 
than ~ kT\,b = 0.2 keV, then the gap voltage is reduced. 

We consider the collision between a gamma ray and a thermal X-ray that is emitted directly from the magnetar surface, 

7+X^e + + e". (37) 

The gamma ray is created with an energy 3 u> 3> 2m e , but below the threshold for single-photon pair creation, to sin 9f, B < 2m e . At 
the point of pair creation, it is nearly completely polarized in the ordinary mode, due to the rapid splitting of extraordinary mode 
gamma rays. The target photon has an energy ux *C m e . The kinematic threshold for pair creation through the channel (|37T i is 
obtained from the equations of conservation of energy and momentum parallel to B, 

w+Ux=(Pt+ m 2 ) 1/2 + (P 2 + m 2 e + 2neB) 1/2 ; 
^uj + ^ x uj x = P + + P-, (38) 

where /i = cos 9ub, = cos 6kB,x are the direction cosines of the gamma ray and X-ray, and P± are the momenta of the final 
electron and positron. One finds 

2tMJx(l-nnx) > (m, + £ „) 2 -w 2 (l-A (39) 

where Eq„ = (m 2 + 2neB) x I 1 and w(l - ^l 2 ) x I 2 = lj sin Bm is the perpendicular energy of the gamma ray. 
Two particular cases are of interest here: 

1. — The gamma-ray remains well below the threshold for pair creation, fx 2 ) 1 ' 2 <C 2m e , so that vacuum polarization effects 
introduce very little dispersion into the relation between u> and k. The target photon is tightly collimated about B in the frame 
where the net momentum of the two photons along B vanishes (0' kBX ~ 2m e c 2 /hu>; Fig. [T). Pair creation with the both final 
particles in Landau state n = has the same threshold condition 

c^ x (l-^) = 2m 2 , (40) 

where we have set /.i ~ 1 . This is identical to the threshold condition in the absence of a magnetic field. 

We have evaluated the cross section for this process to the lowest order in perturbation theory. The tree-level Feynman diagrams 
are calculated using wavefunctions for the electron and positron that are good spin eigenstates (Melrose & Parle 1983a,b). Details 
are provided in the Appendix. We find 



V±\ 3(i T ( B \ (E± 



-2 



2 



(41) 



(E±+P±) 2 + m 2 e + 2eB (E±-P±) 2 + m 2 e +2eB ' 

where E± and P±= P_ = -P+ are energy and momenta of the outgoing electron and positron (both in the lowest Landau state) 
in the center-of-momentum frame, and V± = \P±\/E±. The polarization vectors of the two colliding photons in this frame are 
denoted by = e x k ± ie\. Expression ( HTb simplifies to 

~) 0*7 = ^ 77^2 m-^f (V ± « c) (42) 



c J ' 16 (1+B/B Q ) 

In this section, we employ natural units in which Ti 
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FIG. 1 . — Collision of a gamma ray with a target X-ray. We consider two cases: i) the collision occurs close to the point of emission of the gamma ray, so that 
the gamma ray is well below the threshold for single photon pair creation (fujsmOkB *C 2m e c 2 ); and ii) the gamma ray is just below threshold (eq. 1431 ). In case 
i) the two photons move nearly head on in the center-of-momentum frame parallel to B. In case ii), the cross section is strongly enhanced if the target photon has 
an energy ~ (B/B Q )m e c 2 in the frame where the gamma ray moves perpendicular to B. 



near threshold. The coefficient works out to 0.03OY in a magnetic field B = 4Bq. This result differs somewhat in the dependence 
on B and the numerical coefficients from those of Kozlenkov & Mitrofanov (1986), who used the wavefunctions for the Dirac 
particles that are not good spin eigenstates and did not explicitly evaluate the polarization dependence. 

2. — The gamma-ray is just below the threshold for pair creation, 

cj(1-^ 2 ) i/2 = 2m e -AE. (43) 

In the case of adiabatic conversion of a gamma ray to a pair, AE measures the amount by which the dispersion curve of the 
photon bends below the rest energy of an electron-positron pair. For example, AE ~ 15 keV when the perpendicular momentum 
of the photon is k± = 2m e in a magnetic field B = 4-Bq (eq. 11271 ). 

In this situation, there is a very low threshold energy for the creation of a pair with both particles in Landau level n = 0. 
Transforming to the frame 4 in which the gamma ray moves perpendicular to B, one has 

uo' x > AE. (44) 

This energy is even smaller in the stellar frame, uj x = ^'x/li^ ~f J >x)- When one of the final state particles is in a higher Landau 
state n, one obtains 

uj' x = AE+-(n — + E 0n -m e ), (45) 
2 \ m e J 

which reduces to AE + eB/m e in a sub-QED magnetic field. 

The cross section is strongly enhanced if the virtual electron that mediates the collision between the two photons is on mass 
shell. This is possible only if the energy of the target X-ray satisfies a resonance condition. (The main difference with the 
resonance encountered in cyclotron scattering is that the virtual electron can remain in the lowest Landau state.) We take the 
gamma ray (parallel momentum //u/ = 0) to share a vertex with the outgoing positron (see Fig. @). The parallel momentum Pj of 
the virtual electron is related to that of the target X-ray (p! x oj' x ~ uj' x ) and final-state electron by 

P I =-u' x +P'_ = -P' + . (46) 

The resonance conditions are 5 

{u x -E'_f = (u'-E' + f=P}+ml. (47) 

4 Obtained by a boost with Lorentz factor 7=1/ sin^g ~ 2m e /AE along B. 

5 We can equally well exchange the electron and positron in the final state; both Feynman diagrams contribute to the resonant cross section. 
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We are searching for the lowest-energy resonance, and so take the virtual electron to be in the lowest Landau level. The condition 
of energy conservation 

lu' + lu' x =E' + + E'_ (48) 
is satisfied only if one final particle is in the lowest Landau state, e.g. n+ = 0. Then the solution to equations ( f46b . ( |47| > is 

E - = ^{ 1 + 2^) +n -£>, (49) 

where we have approximated u>' ~ 2m e . Applying eq. d48l shows that the energy of the target X-ray must satisfy the same 
resonance condition as for absorption by a real charge at its first Landau transition, 

, eB 

u' x = -f(l-ii x )uJx = n- . («_>0). (50) 

m e c 

The resonant target photon energy is essentially the same as the threshold energy (|45l l when B < Bq, but lies substantially above 
threshold in super-QED magnetic fields. 
It should be emphasized that this resonance is only approximate when AE > 0. Substituting n+ = into eq. j49l , one sees that 

E' + = m e - l -AE. (51) 

A true resonance is not possible, because one particle in the final state would not be able to satisfy the mass shell condition. But 
an approximate resonance is present when AE <C 2m e . The cross section derived in the Appendix is 

- = \e-^l«(^-\\e?f. (52) 
a j 8 \ur - 2E + / 

In this expression, the frequency a/ and polarization e z of the incident gamma ray, and the energy E' + of the outgoing positron, 
are evaluated in the frame where the gamma ray perpendicular to B. 

The momentum of the positron can be evaluated in terms of the displacement of the energies of the gamma ray and target 
photon from the resonant values u' = 2m e and u>' x = eB/m e : 

i e -oj ) + uo x . (53) 

' m e 

The cross section peaks at P' + = 0, where the denominator of eq. d52l has the minimum value AE 2 . 

The cross section exhibits additional localized peaks where the energy threshold for excitation to higher Landau levels is just 
satisfied (Kozlenkov & Mitrofanov 1986). We ignore this effect in this paper, being interested in the case where the fluxes of 
gamma rays and target X-rays are both cut off exponentially at high frequencies. 

4. HIGH PAIR DENSITIES IN A TURBULENT FLUX ROPE 

A twisted magnetosphere is subject to current-driven instabilities, which can drive fluctuations in the twist and facilitate its 
redistribution across magnetic flux surfaces. Before considering the time evolution of the twist in ^5] we first consider the 
influence that these twist fluctuations would have on the pair density in the outer magnetosphere. We describe a novel state 
for the open-field circuit with very strong and efficient particle heating. Pair creation is driven by low-energy particles and is 
concentrated close to the star, although the heating is distributed over a wide range of radii (see Fig. |2). 

Analogous effects are believed to occur in the magnetosphere of one star in the double pulsar system J0737-3039, which is 
exposed to the relativistic wind of the second, faster pulsar. The dynamic twist can, in both cases, have important consequences 
for the particle energy distribution in the circuit (Thompson & Blaes 1998; Lyutikov & Thompson 2005). On some flux tubes, 
the current that is drawn from the surface of the star will intermittently be pushed above J = pcic, which has the effect of igniting 
pair creation in order to supply the charges that are demanded by the circuit. 

Yet higher current densities will result from the non-linear coupling between torsional Alfven waves propagating in opposing 
directions along the magnetic field, which can create a high-frequency spectrum of modes. The necessary condition for this 
cascade to occur is that the magnetic field lines are sheared over a scale (transverse to the magnetic field) that is somewhat 
smaller than the radius. The cascade is not restricted to the closed magnetic field lines, since an outgoing wave is supplied by 
reflection from the surface of the star. 

Consider a torsional wave of a moderate amplitude 

SB^-B^ (54) 
which propagates along the dipolar magnetic field B = Bns('"/^ns)~ 3 - Here 

B^= ( — ) (—) ' B(r) (55) 









K) 







S PGJC/ \ C 

is the static toroidal field at an angle 8 ~ 6> open from the magnetic axis, where 



y open 



<')=(^j (56) 
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FIG. 2. — The excitation of a torsional wave on the extended magnetic field lines (open and closed) heats the enclosed electrons and positrons when the 
wave energy cascades to a high wavenumber. The fluctuating component of the current diverges in the cascade. There is a critical wavenumber above which 
the relativistic plasma particles can no longer conduct charge at the rate demanded by the waves, and the wave energy is transferred to the particles. When the 
particle density is much higher than jiqj and the plasma is nearly charge-balanced, then waves are excited close to the star, where resonantly scattered photons 
can convert directly to pairs off the magnetic field. The loss of charges from the ope n magnetic field lines can be balanced by pair creation within the innermost 
part of the circuit if the spin period is shorter than the critical value given by eq. i79i . 



is the polar angle subtended by the open bundle of magnetic field lines. In addition, / is the current density that is drawn from the 
surface of the star along the open field lines, pqj = -fi ■ B/2irc is the corotation charge density (Goldreich & Julian 1969), and 
is the spin vector of the star. The wave is decomposed into ingoing and outgoing components, 

(r,0 = SB-Jr+V w t) + SB + Jr-V w t) (57) 



which move with a phase speed V w close to the speed of light. We are interested in the case where the charged particles that 
support the wave have a one-dimensional relativistic distribution, which extends up to some Lorentz factor 7 max . Then 

- /' 1+ lV / 7 1+ Ji£LV'^ 



c k\\c V a 7 V 47 max cjp/ 

In this expression, ujp = (47r«±e 2 /m c ) 1 / 2 is the non-relativistic plasma frequency, and a ~ B 2 / 8n"/ mRX n±m e c 2 . 

The coupling parameter that determines the strength of the interaction between two colliding Alfven waves is k±^ = 
k±8B$/k\\B, where £^ is the rotational displacement of the magnetic field lines and k± and k\\ are the components of the wavevec- 
tor perpendicular to and parallel to B. A cascade develops if the shear of the field lines at the forcing scale is stronger than 

^o^f^V 1 . (59) 



V B 

Here SB^q is the r.m.s toroidal field at the forcing scale, which is take to be 

2tt cP 

— =L~—. (60) 

Km Z7T 

This coupling parameter has been hypothesized to be regulated to kx^ ~ 1 throughout the inertial range of the wave spectrum 
(Goldreich & Sridhar 1995). Then a wavepacket undergoes a significant distortion during a single collision. If the cascade is 
characterized by a constant flux of energy toward higher wavenumbers, 

~ — -k»c = const, (61) 

dt Sir " 

3/2 

then the Alfven wavepackets become increasingly elongated within the inertial range, k± oc kJ . 
The fluctuation in the current density associated with the Alfven waves increases in magnitude toward higher wavenumbers, 

6J= -Z-kxSBttxk 2 / 3 , (62) 

4"7T - L 

whereas the energy density in the toroidal field fluctuation decreases, 



SB\ 5B\ Q fk\\L 



8tt 8tt V 2tt 



(63) 
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*||,max~ 7; ■ (65) 



The current fluctuations at the damping scale are therefore much larger than the mean current density. When the plasma is cold, 
the damping scale is obtained by balancing the r.m.s. current density of the wave with the maximum conduction current that the 
plasma can supply, 

k^SBij, ~ 4nen±. (64) 

The inner scale of the turbulent spectrum therefore sits at a wavenumber hi max ~ 4iren±/B. 

We are, however, interested in the case where the plasma is dilute enough that the charges become relativistic. The heated 
charges are assumed to move only in the direction parallel to B. This assumption is well justified close to the star, where we 
focus our attention, due to the rapidity of synchrotron cooling. In this case, the particles that support the fluctuating component 
of the current must reverse direction: the current fluctuation due to charges that do not reverse direction is negligible, because 
their velocity undergoes only a small change, Sv ~ c/7 2 . We deduce that the current fluctuation is supported by only a fraction 
of the heated charges, with a net density Sn± < n±, and a characterstic Lorentz factor 7d amp that is smaller than the mean Lorentz 
factor (7) ~ SB^ /8irn±m e c 2 . The associated damping scale is given by 

A-Ke§n± 

It turns out to be possible to deduce Sn± and the damping scale without knowing the distribution function of the charges in 
detail. At the damping scale, the energy A£damp deposited per charge in one wave period is comparable to the kinetic energy 
7damp'w«c 2 of the charges that revese direction, 

A£ damp ~ — ~ Tdamp^e . (66) 

on± at K|| max c 

Here 

^i~fi^M. (67) 

dt 57T 

(Except for the primary Goldreich- Julian charge flow, higher-energy particles come in pairs of opposite signs that do not gain or 
lose a significant amount of energy.) Combining eqs. d65l l - d67| ) then gives the damping scale 

2tt /&Tr 2 eB NS R NS \ 1/2 /flR NS \ 2 

fc||,max ~ S~ — > (6«) 

where Z?ns is the polar magnetic field strength. The values of 7d am p and feii max can be determined independently by the physics of 
pair creation, which leads to a constraint on the spin period ( 34. 11 1. 

A comment about Landau damping is in order here. We are interested in the regime where the fluctuating component of the 
twist is comparable to the mean twist, so that a ~ 1 . The Alfven waves then resonate only with the low-energy tail of the particle 
distribution, and the Landau damping time is much longer than the wave period. Since the lifetime of a wavepacket is comparable 
to the wave period in a critically balanced cascade, the effects of Landau damping can be neglected. As a result, Landau damping 
of the Alfven waves is ineffective and the inner scale of a turbulent spectrum is located at the scale where the waves become 
charge starved. 6 

4.1. Maximum Spin Period for High Rates of Pair Creation 

We now consider the range of spin periods and dipolar magnetic fields that can support rapid pair creation with the circuit 
model just described. We focus on the resonant scattering of thermal X-rays by relativistic charges, which requires only modest 
particle energies and can therefore be effective at high particle densities. The creation of pair-creating gamma rays by this process 
is, however, only effective within a relatively small distance from the star. 

The rate of energy deposition by a turbulent cascade is conveniently normalized to the corotation particle density «gj = 
VLB/lirec, so that (1 /nciidU^/ 'dt is approximately independent of radius. If pair creation is concentrated close to the star, 
then n± oc «gj °c r~ 3 on the open magnetic field lines, and the heating rate per electron or positron will be approximately constant. 
Heating is cut off inside a minimum radius 

R«m~-j^—<iL (69) 

^||, max 

where the energy density in the fluctuating toroidal field becomes small, SB^/Stt oc r 2 . 

An upper bound on j? m j n is obtained by requiring that the magnetic field be strong enough for pair creation. One therefore 
requires a short damping scale and a high pair multiplicity. If the spin frequency of the star is too low, then the damping energy 
A£damp becomes too small for the heated particles to produce pair-creating gamma rays. The energy threshold for pair creation 
is lower for ingoing charges, for which the resonance condition (0 is most easily satisfied. We therefore focus on pair creation 
by photons that are resonantly scattered toward the star. 

Pairs that are created by this process can be supplied to the outer magnetosphere only if they reverse direction, that is, only if 
their kinetic energy is smaller than A£damp = 7damp'«,?c 2 (eq. l66l ). The energy of a resonantly scattered photon can be related to 
the kinetic energy j les m e c 2 of the scattering charge via 

„2 ^ ,( B \ „ 2 



£7 = /recoil 7res'«eC ~ ^ — J -y Ks m e C . (70) 
' See Thompson (2006) for a description of the conditions in which these two damping regimes will be found in a cold plasma. 
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(It will turn out that pair creation is concentrated where B < Bq, and so we approximate /recoil ~B/Bq.) The energy of the gamma 
ray is approximately divided in two by pair creation, which leads to the relation 

B{R m in) , n , . 

7damp = ~TTZ 7res ■ ( ' 1 ) 

IBq 

Pair creation in the zone r ~ R m i n will balance the loss of charges from the circuit if 7 res (and therefore the energy AEdamp) 
exceeds a critical value. The resonant scattering rate ( TT7| > depends exponentially on 7 res when the target photon is drawn from the 
high-energy tail of the black body distribution. The value of 7 res that yields a significant rate of resonant scattering decreases as 
one moves away from the star, but its value at a fixed radius depends only logarithmically on the pair creation rate through the 
parameter 2< (eq. ffT8l ) 

B(R min )/B Q 

7res = ~ q ■ (72) 

ZW<jUbb 

One requires ui^ ~ 12 if «±/«gj is almost uniform throughout the open-field circuit. 

The values of R m i n and 7d a mp can now be determined by requiring that photons of energy ( F70b able to convert to pairs off the 
magnetic field. The corresponding condition is 

E 7 r f n > 2m e c 2 , (73) 



where 



Rc(r) / r , 

-jr^U-i (74i 



is the radius of curvature of the field lines, and is assumed to be dominated by the quadrupole/octopole component of the magnetic 
field [see eq. (A9) of Thompson 2008]. One finds, 

^min -4/8 1/4 (^<\ V% ( k B T bb 

— <3.2 £3 Vis^ ) (.O^kevJ ' (?5) 

and 

B(r min) > 3.1 x 10V/»< U (^f) 3/8 (o|^) V8 G - < 76 > 
The corresponding energy of the pair-creating charges is 



>30^ 3/ Xi 4 15 (^l \-^-\ , (77) 



~ \ -5/8 / , „ v 



12 7 \0.5kevi 



>~ \ -1/4 / , t x -1/4 
,1/2 —3/4 I W< \ / KB^bb \ 



and the energy of the resultant pairs is 

7da m p(« ml n)>10B;' s : i5 ^^^j ^^f^vj ■ ( ?8 > 

Charges can be supplied self-consistently to the circuit only if the spin period is shorter than a critical value, which is obtained 
by substituting expressions d75l ) and ( fTBl into eq. (f69t and then into eq. 



P<P crit = 7.7ef < 4 15 < s 4 6 s. (79) 

When the spin period is larger than d79b . the cascade deposits energy too far from the star to deflect ingoing pairs from their 
creation zone at r ~ R m i n back out into the circuit. This expression has the attractive feature of depending only weakly on the 
polar magnetic field Z?ns an d tne degree of curvature of the magnetic field lines; it is also independent of Tbb and uj < . However, the 
normalization does depend on the numerical coefficient (eq. l69l ) relating the minimum radius for pair creation to the damping 
scale: one finds P cr j t oc (£|| jm ax^mm)~ 1,/2 - The electric potential fluctuations at the damping scale will, in practice, have some 
distribution about A£damp/ e ( ec l- Il66l ). with an exponential cutoff at large values. As a result, the supply of pairs to the outer 
circuit is possible if R m [ n is somewhat smaller than the Alfven wavelength at the damping scale, as is assumed here. 

The charges flowing outward toward the open end of the magnetic flux tube continue to gain energy, because the cascade energy 
is deposited broadly throughout the circuit. The gamma-ray emitting particles at r ~ R m i n have an energy 7 res > 7damp (eq. i PTD ). 
and therefore have typically been heated at a larger radius and then partially deflected back toward the star. Without prescribing 
in detail the pair distribution function, it is possible to set a lower bound on the multiplicity in the circuit, 

, , 2iten+c 2?te 5n+c 

m ± = ^iT > ^b^- (80) 

Substituting eqs. (l65l l and ( |69l , one finds 



1 f flR min \ 1 1n 3—l/8 n -l/4 ( P\ ( ^B^bb V^ 8 ,,,, 

M±> 2[—) =4xl0£3 ^"(elJlastev' ' (81) 
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Given a constant outward flux of pairs with a radial velocity comparable to c, the mean energy per charge is 

r dU cas 



(7(0) ■ 



M±n G s dt 

P V 3 / k B T bh V 1/8 fSB^oV ( J 2 



-|4~l/8 R 5/4 3 

fc 3 b ns,i5"n; 

Note that (7) > 7d am p at r > R m i n : most of the heating is at the low-momentum end of the distribution function f(P). If df/dP ~ 

—1/2 

at low momentum, then 6n±/n± ~ 7dam P / (7)- But Sn± oc 7 da 4 p , and so we deduce 

\ 2/3 

r 



7damp(0 = 7damp(flmin) ^nR ) ' ^ 

The damping scale 27r/A:|i „ oc 5ri± increases as the i power of the radius. 

A competing circuit solution that could lead to such large pair multiplicities is an outer gap bounded by the surface fi • B = 
(Cheng, Ho, & Ruderman 1986; Cheng & Zhang 2001). The voltage across such a gap can approach the total open-field voltage. 
Seed charges can be created in the gap by collisions between curvature gamma rays and thermal X-rays from the neutron star 
surface. The energetic charges created in such a gap will, however, be beamed down toward the star, and pair creation will be 
concentrated in a zone where the open field lines are occulted by the star. Such a voltage structure is therefore not a promising 
source for the observed emission of the radio magnetars. 

5. VARIABILITY IN THE OPEN-FIELD VOLTAGE DUE TO MAGNETOSPHERIC ACTIVITY 

The variations in the spindown torque of SGRs and AXPs are remarkable in their amplitude and longevity: the spindown rate 
is observed to increase smoothly over a period of months, by up to a factor ~ 4 (Kaspi et al. 2001; Woods et al. 2002; Woods et 
al. 2007; Tam et al. 2007; Camilo et al. 2007a, 2008). The strength of the torque variations correlates roughly with the level of 
activity that the magnetar has sustained over the preceding years, but this correlation is only indirect. In several cases the torque 
increase is observed with a delay of months following an episode of X-ray burst activity, and can persist long after the intial active 
period. The burst emission of SGRs and AXPs is highly intermittent (Woods & Thompson 2006): a nice example is provided 
by the AXP IE 2259+586, whose activity as a hard-spectrum X-ray source was localized within a day, and was followed by an 
extended cooling of the thermal X-ray emission (Woods et al. 2004). In the most active sources, the magnetar can remain in a 
new spindown state for years at a time, with the torque several times larger than in the preceding state. 

Strong variations in spindown torque, extending over a period of at least several months, have also been detected in the radio 
monitoring of the transient AXP XTE J1810-197 (Camilo et al. 2007b). This torque variability persisted long after the initial X- 
ray activation (in fact, after the hard-spectrum component of the transient X-ray emission had largely decayed away). Even more 
remarkably, the torque increased by at least a factor ~ 3 compared with older X-ray timing data, making the relative increase in 
the spin down rate comparable to that seen in the most active SGRs. 

These torque variations have been interpreted in the magnetar model in terms of a restructuring of the magnetosphere. Two 
effects have been considered which can increase the open magnetic flux, the open-field current, and the torque acting on the star. 

1 . - A persistent static twist may be sustained in the closed magnetosphere, driven by the unwinding of an internal toroidal 
magnetic field (Thompson et al. 2002; Beloborodov & Thompson 2007). In the process, magnetic helicity is expelled from the 
interior of the star. The twist is stabilized as long as the outward flux of helicity across the outer boundary of the magnetosphere 
is kept small. 

As a dipolar magnetic field is twisted, the poloidal field lines flare out slightly. The large size of the corotating magnetosphere 
then allows a significant increase in the open-field current. In a first approximation, one additional parameter is introduced into 
the description of the magnetosphere, namely the twist angle A<f> on the field lines that are anchored near the magnetic poles. In 
the self-similar construction of Thompson et al. (2002), the radial index of the magnetic field softens from B oc r" 3 to B oc r~ K 
with K < 3. The angular width of the bundle of open field lines therefore increases from 9l p( . n = Qr/cX.0 1 

Cn=(^)" 2 (*<3). (84) 

The voltage across the open flux bundle is proportional to 6^ pen . A twist of one radian forces a reduction in K to 2.85 (see Fig. 2 
of Thompson et al. 2002). At a fixed rate of spin, the voltage increases by a factor (cP /2nR NS ) 3 ~ K = 4.7 (P/6 s)° A5 R^6' as does 
the polar magnetic field that is inferred from the magnetic dipole model. 



2. - A persistent relativistic wind is driven by internal seismic activity within the star. This wind combs out the dipolar field lines 
into a radial configuration in the outer magnetosphere (Thompson & Blaes 1998; Harding et al. 1999; Thompson et al. 2000). 

The conversion of internal magnetic and seismic energy to a persistent wind has an uncertain efficiency: it involves the transfer 
of energy across magnetic flux surfaces from an injection zone situated far inside the speed-of-light cylinder. For example, 
torsional waves in the outer magnetosphere have a much lower frequency than do internal elastic modes of the neutron star crust: 
the fundamental toroidal mode of the crust has a period P e iastic ~ 0.02 s (Strohmayer et al. 1991) and couples to dipolar field lines 



7 This expression is accurate when K is only slightly smaller than 3; it neglects any change in the distribution of flux across the neutron star surface. 



14 



THOMPSON 



with a maximum radius R mdx ~ 0.4cP e i ast i c , only ~ 10~ 2 (P/6 s) of the light cylinder radius. A wind can therefore be driven by 
the excitation of waves in the closed magnetosphere, but the power dissipated on field lines that open out across the light cylinder 
may be only a small fraction of the total. 

Measurements of the torque behavior of magnetars offer ways of discriminating between these two models. 

A transient twisting up of the magnetic field has the effect of increasing the open-field voltage and reducing the critical spin 
period below expression (1791 . If the detection of the radio emission from an individual magnetar typically requires a broad radio 
beam, and therefore depends on a high pair density on the open field lines (see Thompson 2008), then one can understand why 
an individual source would become visible as a radio pulsar following a period of X-ray burst activity. Starting with a value of 
f cr it larger than the present spin period, a twisting of the magnetic field can force P cr [ t to a smaller value and ignite a high rate of 
pair creation on the open field lines. 

The delays that are observed between burst activity and torque increases can be interpreted in terms of the resistive evolution of 
the current within the closed magnetosphere, which has a timescale of ~ 1 yr close to the star (Beloborodov & Thompson 2007). 
Their existence is less clear in the case where the torque change is driven by persistent seismic activity. An attractive feature 
of the twisted dipole model is that small injections of twist from the neutron star interior can force large structural changes in 
the outer magnetosphere (where the magnetic field energy is relatively small). Indeed, the energetic output of XTE J1810— 197 
during its period of X-ray activity was modest by the standards of the Soft Gamma Repeaters; but the amplitude of the subsequent 
change in spindown torque was as large as is seen in the most active SGRs. 

The measured upper bound on the spin periods of magnetars places additional strong constraints on how the magnetophere is 
restructured. A persistent wind of particles and Alfven waves drives an exponential increase in the spin period, 



dP 4^i NS / 

^wind 

dt 3/ N s 



(85) 



where /i N s = ^Bns^ns * s tne ma g ne tic moment. 8 This formula applies when the wind luminosity L w j n d is much larger than the 
magnetic dipole luminosity Lmdr = §Mns^ 4 / c3 - One finds 

^MUd =60 (P\ 2 m 



dP/dt\ M DR #NS,15^NS,6 \6s 

We deduce that a relativistic wind can dominate the spindown torque of a magnetar only rarely — otherwise the spin periods 
of magnetars should cover a broad range. The efficiency of this torque mechanism does not depend in an obvious way on the spin 
period and the rate of pair formation the open magnetic field lines. Alfven waves in the closed magnetosphere are supported by 
charges particles of both signs, that are drawn directly from a light-element surface layer, or supplied by pair creation near the 
magnetar surface (as in the circuit model of Beloborodov & Thompson 2007). 

5.1. Timescale for Resistive Evolution of the Closed-Field Current and the Open-Field Voltage 

Strong variations in radio brightness and pulse shape are detected from the two radio magnetars on timescales of minutes to 
hours (Camilo et al. 2007b, 2008). These variations have a possible interpretation in terms of the resistive evolution of currents 
in the outer closed magnetosphere, which drives variations in the open-field voltage. 

Let us first consider the gradual redistribution of twist that would result from smooth gradients in the field-aligned voltage 
across the poloidal flux surfaces. To estimate the timescale in the outer magnetosphere, we consider the power dissipated on field 
lines that extend to a certain maximum radius R mdx . The self-similar model of a twisted magnetosphere described in Thompson 
et al. (2002) has a current density of the form 

J = aB; a= — sin 2 f? . (87) 

47T/? NS 

Here 6*o is the polar angle at which a given field line intersects the star, and Acf> is the net twist angle between the two magnetic 
footpoints. We assume that A<j) is constant in the outer magnetosphere, although in general it will vary across the poloidal flux 
surfaces. Then the toroidal magnetic magnetic field is given by 

A6sm 3 9 

B^r,6)= ^ 3 MNS , (88) 
and the net current flowing across a radius r > /?ns is 

^ 2 / r V 2 

— cBkkRxsk 



/(> r) = j J,.(r, 9)2irr £ smOdd = cB NS R m ■ (89) 

Given that each closed flux surface develops the same voltage the dissipation rate is 

L d U>r) = I(>r)$ e . (90) 

8 In the case of a centered dipole, /xns can be related to the polar magnetic field Bns and the average surface field (#ns) via fifis = ^5nsRJj S = (Bns)^S' 
It is (Bims) that is usually quoted in the pulsar literature. Hence a polar magnetic field B^s = 10 15 G, the usual normalization in this paper, corresponds to 
(Bns) =5 X 10 14 G and a magnetic moment /jns = 5 x 10 32 S^ s 6 G cm 3 . 
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The decay time of the twist can be deduced from Ldiss and the energy in the toroidal field, which is 

E 4 >r )s f (£f . (9D 

The timescale for the resistive evolution of the twist at radius r is 

Here we have normalized r to the light cylinder radius R\ c = c/fi. The field-aligned voltage <£> is the integral of - En = —E-B along 
a closed magnetic field line between its two footpoints, 9 The value of $ in the outer magnetosphere depends on details of the 
return current toward the star: a gap forming at the surface of an X-ray bright magnetar has a voltage e$ ~ 1 GeV (Thompson 
2008), which is similar to the voltage deduced from the 1-dimensional circuit model of Beloborodov & Thompson (2007). At 
constant $, the resistive timescale increases toward the star. 

Our evaluation of f res is 1-2 orders of magnitude longer than the shortest timescale for radio flux variability observed in the 
magnetars XTE J1810-197 and IE 1547.0-5408. The underlying voltage fluctuations may therefore be triggered by current- 
driven instabilities in the outer magnetosphere. In the next section, we consider the origin of these instabilities, and how they are 
regulated by pair creation. 

5.2. Stability of the Magneto spheric Twist 

A substantial increase in the open-field current and spindown torque will result from a twisting of the magnetosphere only if 
the twist is distributed over a broad range of radii. If the twist were concentrated within a distance ~ 2R^s of the neutron star, the 
effect would be much more limited. The detection of strong nonthermal X-ray emission from a magnetar does not necessarily 
imply that the source is in its highest torque state. For example, the AXP 4U 0142+61 is a copious source of 100 keV X-rays 
(Kuiper et al. 2006), but its characteristic age is relatively long (? sd - = P/2P = 70 kyr); its outer magnetosphere may therefore be 
only weakly twisted. 

How stable is a magnetosphere with a broadly distributed twist? A self-similar configuration with a constant value of the twist 
angle is not in its minimum-energy ('Taylor-relaxed') state. The constant of proportionality a between the current density J and 
magnetic flux density B varies with position (eq. l87l ). In particular, it depends on the polar angle 9q at which a given field line 
intersects the star, which provides a convenient label for the poloidal flux surfaces. In what follows, we allow A<f> to be a function 
of O . 

Such a configuration with a decreasing toward the magnetic symmetry axis is susceptible to redistribution of twist into the 
outer magnetosphere. The kinetic energy of the charges that supply the closed-field current is small compared with the magnetic 
energy, which means that the magnetosphere is nearly force free. The energy of a force-free magnetic field is minimized at 
constant a (Woltjer 1958). When A(j> is constant, most of the magnetic helicity is stored in the inner magnetosphere: 

H(r) ~ [irr 2 B P (r)] 2 A<t>(r) oc -^2. (93) 

The same is true when a is constant: in this case A(j> oc 1/ sin 2 6'o oc r, and ri oc r~ l . The toroidal field energy is therefore 
minimized when most of the helicity is concentrated close to the neutron star. At the same time, the redistribution of a small 
portion of the helicity into the outer magnetosphere can drive the twist angle to large values. 

To describe the stability properties of the magnetosphere in more detail, it is useful to define appropriate angular coordinates. 
We focus on axisymmetric equilibria, for which the magnetic field can be written as 

B = Bp+Bt = Vtfpx + Vf^x V(^). (94) 

Here typ = /ZNssin 2 #/r is tne A ux coordinate for a dipolar magnetic field, and = </>/rsin# in spherical coordinates (r,9,<f>). 
The toroidal magnetic field corresponding to eq. (l87l i is given by eq. (l88l . The toroidal flux threading a poloidal flux surface of 
cross-section S that is anchored at angle 9q is 



* (0 O ) = J rdrd6B 4 (r,9). (95) 
The surface integral d95l l can be transformed into an integral over 9 and 9q using the relation r = R^s sin 2 9/ sin 2 9q, 

'o)= / sm9' co&9' Q d9 / sm9d9. (96) 



Rm 

The distribution of toroidal flux with surface polar angle is 

_ 2A0^ NS 
d sin 9q /?ns 



sin 6*0 cos 9q, (97) 



9 This quantity represents the line integral of (l/c)dAu jdt - not the change in the scalar potential, which must vanish given the high electrical conductivity of 
the neutron star. 
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FIG. 3. — Cross section of the current density J and safety factor q in a plane perpendicular to the long axis of a tokamak torus (panels a,b) and to the symmetry 
axis of a twisted dipole magnetic field (panel c). In a tokamak, q must have a minimum near the axis if / does. A given mode then can interact with two rational 
surfaces, which accentuates the growth rate of current sheets at the rational surfaces. The net effect is a rapid flattening of the current and safety factor profiles. 
By contrast, the q profile can remain quite flat near the symmetry axis in a neutron star magnetosphere, even through J has a minimum on the axis. 



and the dual angular variable is 

2tt 2 \ co&8() / 

The strength of the winding in the magnetosphere can be described in terms of the relative periodicities in the <fi and t] directions, 



77 1 / co&9 , 

' 1 7T • (98) 



1 _d<t> (dr)\ l A0(6» o ) 



q(9o) drj BgsinO \d0 J 2tt 



cos0 o . (99) 



A twisted magnetic field is subject to a variety of current-driven instabilities. Ideal MHD instabilities set in at strong values 
of the twist: the kink instability is present when the 'safety factor' q < 1 in cylindrical geometry (e.g. Wesson 2004). On the 
other hand, the flaring of the dipolar magnetic field lines has a strong effect on the spindown torque even when A(j> ~ 0.5- 1, 
which corresponds to q ~ 5- 10. The hard X-ray emission of magnetars can be powered by even weaker twists in the inner 
magnetosphere: only A0 ~ 0.1 is required (Thompson & Beloborodov 2005, Beloborodov & Thompson 2007), corresponding 
to q ~ 10 2 . 

The redistribution of twist by tearing instabilities within the closed magnetosphere of a neutron star must therefore be con- 
sidered. Resistive kink instabilities are known to trigger a redistribution of current across magnetic flux surfaces in tokamak 
plasmas. They can be excited at lower values of the twist, and are sensitive to the current distribution within the tokamak torus. 10 

Here there is an important distinction between a tokamak plasma and the twisted magnetosphere of a neutron star. The toroidal 
current in a tokamak plasma torus is observed to relax rapidly to a constant-a state if the current is concentrated initially near the 
outer boundary of the torus. As a result, the profiles of the current and the safety factor q both become much flatter in near the 
long axis of the torus (Fig. [3^,b). The outer magnetosphere of a neutron star can, by contrast, sustain a current distribution in 

10 Our notation here is the inverse of that in the tokamak literature: the unstable component of the twist is in the azimuthal direction <j>, whereas in a toroidal 
tokamak plasma it is in the 8 direction. The poloidal magnetic field of a neutron star is supported by currents embedded in the star, and the external toroidal field 
by currents flowing through the magnetospheric plasma. In the case of a tokamak, the toroidal magnetic field is supported by external currents that loop around 
the long axis of the torus, and the poloidal field (Bg) is supported by plasma currents. 



ELECTRODYNAMICS OF MAGNETARS: PAIR CREATION AND DYNAMIC PARTICLE HEATING 



17 



which the g-profile is flat across the poloidal flux surfaces (A0 is constant) even though the current density decreases toward the 
magnetic axis, J oc 9 2 (Fig. [3};). 

This difference in the shape of the ^-profiles in unrelaxed tokamak and neutron star plasmas has interesting implications for the 
excitation of resistive kink instabilities. The tearing modes which facilitate the redistribution of twist across the (toroidal) flux 
surfaces are concentrated on rational surfaces where 

(B .V)<5J=^(^-«^J = 0. (100) 

Here SJ oc e l( - m0 ~"^ is the current perturbation and R is the large radius of the plasma torus. The analogous form for the current 
perturbation in a neutron star magnetosphere is 



SJ oc exp 



, , cosfl 
imm—mn 1 



cos 6*o / 

where the roles of the 9 and cj> angular variables are reversed. It is straightforward to check that 



(101) 



(B-V)tfJ= — f — -m-n ) SJ, (102) 



r \2tt 

where we have set cos#o — > 1 for the most extended dipole field lines. 

When the current distribution across the long axis of a tokamak torus has a localized maximum, then so does the distribution 
of q. As a result, a single resonance q = m/n can appear at two separate surfaces within the torus (Fig. [3). These double 
tearing modes appear to be important in the process of Taylor relaxation (e.g. Stix 1976). They can be avoided in the twisted 
magnetosphere of a neutron star if the twist has the same sign on the open and closed magnetic field lines in both magnetic 
hemispheres. The presence of a highly conducting boundary also tends to suppress tearing modes at rational surfaces near the 
outer boundary of tokamak plasmas (e.g. Wesson 2004). 

The spin of the star induces a twist on the open magnetic field lines with the same sign as the twist that is sustained by the 
outward Goldreich-Julian current. A reversal in the sign of the twist on the open field lines, if sustained over many rotation 
periods, therefore demands a reversal in the sign of the conduction current. This is possible only if a polarizable plasma is present 
on the open field lines: the space charge density on the open magnetic field lines just above the surface of the neutron star must 
be the same as pcj- On this basis, Thompson et al. (2002) conjectured that the rate of reconnection near the magnetospheric 
boundary would be sensitive to the pair production rate in the open-field accelerator. 

The stability of the twist in the outer magnetosphere depends on whether the safey factor q (eq. 11991 ) can maintain a flat 
profile across the open-field boundary. In the standard pulsar model, the open-field current in each hemisphere is compensated 
by a return current of the same magnitude that flows along a narrow annulus of closed field lines just inside the magnetospheric 
boundary. If the open-field current had the opposite sign to the Goldreich-Julian current in one hemisphere, then the toroidal field 
would go to zero on the flux surface bounding the return current sheath. This would have the effect of allowing double rational 
surfaces to form, in close analogy with a Tokamak plasma (Fig. [3]). By contrast, when the open-field current is able to reverse 
sign, then the ^-profile can remain flat and double rational surfaces do not form. 

6. CLUSTERING OF SPIN PERIODS IN THE MAGNETAR POPULATION AND 
THE SUDDEN BRIGHTENING OF RADIO MAGNETARS 

X-ray bright magnetars 11 range in spin periods between 0.33 and 11.8 seconds. The upper cutoff to the spin distribution 
is statistically significant (Psaltis & Miller 2002). The distribution is also quite strongly peaked: removing two sources (IE 
1547.0-5408, P = 2s, and PSR J1846-0258, P = 0.33 s), the remaining spin periods lie between 5.2 and 11.8 s. If the stars 
underwent simple magnetic dipole spindown, then a range of 2.3 (6) in spin period would correspond to a range of ~ 5 (~ 30) 
in age. The measured characteristic ages actually have a much wider distribution, from several hundred years in the most active 
SGRs to 2 x 10 s yrs in the AXP IE 2259+586. The observation of strong fluctuations in spindown torque in individual sources 
shows that part of the dispersion in characteristic ages is due to variations in magnetospheric structure. This measured torque 
variability cannot be related a narrowing of the spin period distribution unless the underlying mechanism is sensitive to the rate 
of pair creation in the outer magnetosphere. 

We are interested in the case where the rate of injection of toroidal field energy into the inner magnetosphere (when averaged 
over the episodes of X-ray activity) is much higher than the spindown power. The net twist on the magnetic field lines that 
thread the closed magnetosphere and the interior of the star can be described by the magnetic helicity H = J A ■ BdV, where A 
is the vector potential. The toroidal magnetic field in the closed magnetosphere is susceptible to reconnection with the external 
magnetic field. We explained in $5.2l how this process could be impeded by a high rate of pair creation on the open magnetic 
field lines (see also Thompson et al. 2002). 

A slow loss of helicity across the outer boundary of the closed magnetosphere can force large-amplitude fluctuations in the 
twist on the open magnetic field lines, and reverse its sign in one hemisphere. The required rate of loss is 



dU 
dt 



[irRlB P (R lc )] 2 Q, (103) 



" Observed with persistent luminosities of 10 35 erg s 1 or above over a period of years or decades. See the McGill Magnetar web page for an up-to-date list 
of magnetar properties; http : / / www, physics .mcgill ■ ca/~pulsar /magnet ar/main . html| 
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where R\ c = cP/2tt=3 x 10 (P/6 s) km. The corresponding damping time of the twist A<j) close to the star is therefore 

A6(r) Ad, ( P\ 3 

~:30— - — yr. (104) 



\dA^)(r)/dt\i c R NS ,(, \6s 

In the case of a magnetar, the damping time in the inner magnetosphere is comparable to the resistive decay time of several years 
(eq. 192)). 

We have shown in §|4]that a very high pair density can be sustained in the outer magnetosphere if a dynamic twist is excited and 
then damped through a turbulent cascade to high spatial wavenumbers. This process depends on a high rate of particle heating 
close to the star, where the rate of resonant scattering is high, and the resonantly scattered photons are energetic enough to convert 
to pairs. The critical spin period P cr i t (eq. l79l ) for such a gas to be sustained on the open magnetic field lines is several seconds, 
and depends weakly on the polar magnetic field strength and the curvature of the open magnetic field lines. 

6.1. Metastable Twisted Magnetosphere and Torque Decay 

We are led to conjecture that the outer magnetosphere can enter a metastable state in which the twist is lost across the mag- 
netospheric boundary at a controlled rate, just high enough to excite large fluctuations in the twist over a single rotation of the 
star. This metastable state depends on the presence of a dense pair gas in the outer magnetosphere, which can limit the rate of 
reconnection of the closed magnetosphere with the open-field bundle. The detection of hard X-ray emission at a level Lx — 10 35 
ergs s" 1 from AXPs with little known X-ray burst activity (Kuiper et al. 2004, 2006) leads to the remarkable inference that 
magnetic helicity can be stored in the magnetosphere for several years. Hard X-ray emission of a comparable luminosity has 
been detected from SGR 1900+14 in a relatively quiescent state, nearly a decade after the last major X-ray burst activity (Gtitz 
et al. 2006b). Very bright emission (~ 10 36 ergs s _1 ) was detected from SGR 1806-20 by Integral within two years of its major 
2004 outburst (Mereghetti et al. 2005). 

The observed peak in the distribution of magnetar spins can then be explained if most of the spindown torque is accumulated 
in this metastable, twisted state, and not when the magnetosphere is more nearly dipolar. Once the spin period of the star grows 
beyond the value d79l ), the magnetosphere cannot sustain the inflated, twisted state. Of course, there is some hysteresis in thie 
process, because a high rate of pair creation can be reignited by the process of twisting (©, but this only introduces a spread of 
a factor ~ 2 in the value of P cr j t . We also emphasize that this does not exclude an important contribution of a particle wind to the 
torque when the outer magnetosphere is nearly dipolar (Spitkovsky 2006). 

We must therefore examine whether most of the spindown can realistically occur when the magnetosphere is twisted. The duty 
cycle of the magnetospheric twist can be determined by combining the magnetospheric dissipation rate (which can be measured 
in the hard X-ray band) with net twist of the internal magnetic field. The internal field cannot be directly measured, but the 
energetics of the giant flares of the SGRs point to the presence of internal magnetic fields approaching 10 16 G (e.g. Hurley et al. 
2005). 

We suppose that, on multiple occasions, a certain amount of helicity is injected into the magnetosphere due to the release of 
sub-surface stresses, and is stored there. The decay time of the twist can be deduced from eq. d92b . by substituting explicitly for 
the magnetospheric dissipation rate Ld; ss , 

E<p _ 32 Ldiss^NS 

Wy -z d ~-io5^$r 

= 0.1L dlss ,35^NS.6 f^oTy) y r - (105) 

If the current is concentrated on poloidal field lines that are anchored within an angle 8 < 8 C of the magnetic axis, then it is easy 
to show that Ldiss oc 8 A C and f decay oc 6* 2 . 

The net duration twist of the twisted state in the magnetosphere is determined by the winding angle A0, of the internal magnetic 
field. Given that a toroidal magnetic field is present to a depth AR below the stellar surface, one has A<f>j ~ (A/?//?ns)-B^.//-Bns 
and 

, _ A(j)i _ 4 B+jRssAR 

' ^decay — 



Acj) 105 $ e c 



The spindown rate of the magnetar is accelerated when its magnetosphere is twisted. We define 

(i/)a*=Mmb, (107) 
where vq is the frequency derivative in the absence of twisting. The acceleration factor is quite large, N\, ~ 10, if the twist angle 
is as large as 1 radian (Thompson et al. 2002), but observations suggest Nb < 4. Given that f tw i st is smaller than the total age 12 
?sd,o of the magnetar, the relative frequency change due to spindown in the twisted and untwisted states is 

i^± =Nb ^ (108) 
(Az/) fsd.o 

12 We neglect any departures from a simple magnetic dipole torque model, except for those due to twisting. Thompson et al. (2002) found that the braking 
index n = vv / '(u) 2 is slightly smaller than the dipole value (n = 3) in the twisted state, but this feature is not important for present purposes. 
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Consider an SGR with an observed spindown age / s d,o/M> = 200 yrs, which is the minimal value obtained from the spin 
monitoring of SGR 1806-20 (Woods et al. 2007). Most of the spindown will occur in the twisted state, (A^)a0 > (A^)o, if 

( By \ ( $ e y 1 (AR\ 



> 1.7, (109) 



and if the loss of helicity due to non-axisymmetric instabilities in the magnetosphere can be neglected in comparison with ohmic 
damping in situ. If the voltage is as small as ~ 10 8 V, then twisting of the magnetosphere can continue to enhance the time- 
averaged torque as long as the spindown age is shorter than ~ 10 4 yrs. 

6.2. Torque Decay due to Dipole-Spin Axis Alignment? 

In this picture, the spin evolution of magnetars is driven by departures from a dipole structure in the closed magnetosphere. 
The limiting spin period for magnetars is then not directly tied to the radio death line of pulsars. We should, nonetheless, 
consider seriously the possibility that an active pair accelerator has been entirely quenched on the open magnetic field lines of 
many magnetars. X-ray pulse measurements of magnetars would then provide valuable information about the electrodynamic 
properties of pulsars on both sides of the radio death line. 

This intriguing possibility encounters two basic difficulties. First, pair creation by the resonant scattering of thermal X-rays 
is very effective near the X-ray bright surface of a magnetar, and a pair discharge can be maintained on the open magnetic field 
lines at spin periods significantly longer than the observed range of magnetars (e.g. Medin & Lai 2007). Second, a reduction in 
the open-field particle luminosity would result in a large reduction in torque only if the magnetic dipole and rotation axes were 
nearly aligned, so that the magnetic dipole luminosity in vacuum is much smaller than the particle luminosity in the active state 
(e.g. Contopoulos & Spitkovsky 2006). 

Observations of PSR B 193 1+24, which undergoes radio nulls of a month-long duration, show the torque is reduced by a factor 
~ | when the radio emission is off (Rea et al. 2006). If fi and SI are inclined by an angle X1931 in this object, then then one 
deduces that the ratio of magnetic dipole luminosity to particle luminosity is, more generally, 

J^2L^ 2 _^!a_ (110) 

^particle sin X1931 

in a radio-bright state. The characteristic age of the AXP IE 2259+586 is ~ 2 x 10 5 yrs, some 20 times longer than the estimated 
age of the surrounding supernova remnant CTB 109 (Sasaki et al. 2004). To explain a reduction of a factor ~ 20 in torque 
following the quenching of radio emission, one would require that sin 2 x < sin 2 Xi93i/40. 

There are strong arguments against the alignment of fi and SI in most magnetars. The rotational bulge of a slowly rotating 
magnetar is much smaller in amplitude than the variation in the lengths of the principal axes due to the internal magnetic field. 
When the internal field is predominantly toroidal, the star becomes prolate in shape and its two largest principal moments of 
inertia lie in the plane perpendicular to the magnetic symmetry axis. The star reaches a state of minimum rotational energy (at 
fixed spin angular momentum) when SI lies in this plane (e.g. Cutler 2002). If the external magnetic moment were aligned with 
the symmetry axis of the internal field, then there would be no reason to expect alignment between fi and SI. 

A near alignment between fi and SI would be easier if the magnetic moment were distributed randomly with respect to the 
axis of the internal toroidal field. But in that case, the probability of finding alignment to the degree of precision required for IE 
2259+586 is only ±x 2 ~ 0.01. 

7. SUMMARY 

1 . - A gamma ray that is created below the threshold energy for pair conversion can adiabatically convert to bound positronium 
(Wunner et al. 1985; Shabad & Usov 1986). We show that such a positronium atom is rapidly dissociated in the presence of an 
intense flux of thermal X-rays from the magnetar surface, through the resonant absorption of an X-ray by one of the two charges. 
Positronium is also be photodissociated by infrared radiation if the surface luminosity at ~ 0. 1 eV exceeds ~ 10 29 ergs s" 1 . This 
minimum luminosity is only ~ 10" 2 - 10" 3 of the observed infrared output of the AXPs, but some <~ 10 6 times the flux expected 
from the Rayleigh-Jeans tail of the X-ray blackbody. 

2. - We have calculated the cross section for the reaction 7 + 7 — > e + + e~ in the presence of an ultrastrong magnetic field, using 
electron wavefunctions that are eigenstates of the electron spin. We find a resonance in the cross section when the gamma ray 
is, by itself, just below the threshold energy for pair creation. Away from this resonance, the value of the cross section differs 
slightly from previous results. 

3. - We have considered the generation of e ± pairs when the magnetic field lines are strongly turbulent, and have found a self- 
consistent state with a very high pair density and a low mean energy per particle. In this situation, a high-frequency spectrum of 
Alfven waves is created via a turbulent cascade. The current density diverges in this cascade. There is a critical wave frequency 
above which the plasma is not able to support the current, and the energy of the waves is transferred to the longitudinal motion of 
the charges. When the pair density is high, particles are heated close to the star, where resonant scattering is effective. There is a 
critical spin period of ~ 5-7 seconds, below which the pair creation rate in this inner zone is high enough to compensate the loss 
of pairs across the light cylinder and into the pulsar wind. This result does not depend on the details of the particle distribution 
function within the circuit. 

4. - The timescale for the ohmic evolution of the twist in the outer magnetosphere is relatively short (<~ 0.1 day). The fast 
variability that is observed in the brightness of the radio magnetars therefore has a plausible explanation in terms of current- 
driven instabilities on the closed magnetic field lines. 
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5. - We have described in some detail a conjecture that the outer magnetosphere can remain strongly twisted only when the pair 
density is very high on the open magnetic field lines. In this possible metastable state, the loss of magnetic helicity across the 
magnetospheric boundary is just high enough to sustain large fluctuations in the twist, and high rates of pair creation. Part of 
the motivation comes from the consideration of the growth of current-driven instabilities in the outer magnetosphere, which are 
strongly enhanced if the sign of the current on the open field lines matches that flowing into the polar regions of the rotationally- 
driven wind. 

6. - The cumulative torque that can be exerted on the spin of a magnetar in the active magnetospheric phase has been considered, 
and it has been shown that it can exceed the cumulative magnetic dipole torque. For this reason, the observed clustering of 
magnetar spins may be connected to a loss of stability of a twisted magnetosphere at a characteristic spin period of several 
seconds 

As is discussed in further detail in Thompson (2008), the dynamic model of the outer magnetosphere investigated here has a 
promising application to the radio magnetars, in particular to their hard spectra, broad pulses, and relatively high luminosities (in 
comparison with the spindown luminosity). The consistency of this pair-rich circuit solution does not, by itself, imply a radio 
death line for magnetars - a weaker pair discharge appears possible at significantly longer spin periods (e.g. Medin & Lai 2007). 
However, radio emission from the inner magnetosphere would be much harder to detect, due to the effects of beaming. 

This work was supported by the NSERC of Canada. I am indebted to Andrei Beloborodov for explaining some of the literature 
on positronium creation to me. His suggestion that a hybrid gamma-ray /electron-positron pair might be able to photodissociate 
led to the calculation in jjA.21 I also thank Cornell University for its hospitality when part of this work was completed. 

REFERENCES 



Adler, S. L. 1971, Annals of Physics, 67, 599 
Baring, M. G., & Harding, A. K. 1998, ApJ, 507, L55 
Baring, M. G., & Harding, A. K. 2001, ApJ, 547, 929 
Baring, M. G, & Harding, A. K. 2007, Ap&SS, 308, 109 
Beloborodov, A. M, & Thompson, C. 2007, ApJ, 657, 967 
Berestetskii, V.B., Lifshitz, E.M., & Pitaevskii, L.P. 1982, Oxford: Pergamon 
Press 

Camilo, E, Ransom, S. M., Halpern, J. P., Reynolds, J., Helfand, D. J., 

Zimmerman, N., & Sarkissian, J. 2006, Nature, 442, 892 
Camilo, E, et al. 2007a, ApJ, 663, 497 

Camilo, E, Ransom, S. M., Halpern, J. P., & Reynolds, J. 2007b, ApJ, 666, 
L93 

Camilo, E, Reynolds, J., Johnston, S., Halpern, J. P., & Ransom, S. M. 2008, 

ArXiv e-prints, 802, arXiv:0802.0494 
Chang, H.-K. 1995, A&A, 301, 456 

Cheng, K. S., Ho, C, & Ruderman, M. 1986, ApJ, 300, 500 
Cheng, K. S., & Zhang, L. 2001, ApJ, 562, 918 
Daugherty, J. K., & Ventura, J. 1978, Phys. Rev. D, 18, 1053 
Daugherty, J. K., & Harding, A. K. 1989, ApJ, 336, 861 
Dermer, C. D. 1990, ApJ, 360, 197 

Durant, M., & van Kerkwijk, M. H. 2006a, ApJ, 650, 1070 
Durant, M., & van Kerkwijk, M. H. 2006b, ApJ, 650, 1082 
Erber, T. 1966, Reviews of Modern Physics, 38, 626 
Fernandez, R., & Thompson, C. 2007, ApJ, 660, 615 
Goldreich, P., & Julian, W. H. 1969, ApJ, 157, 869 
Goldreich, P., & Sridhar, S. 1995, ApJ, 438, 763 

Gonthier, P. L., Harding, A. K., Baring, M. G., Costello, R. M., & Mercer, 

C. L. 2000, ApJ, 540, 907 
Gotthelf, E. V., & Halpern, J. P. 2005, ApJ, 632, 1075 
Halpern, J. P., & Gotthelf, E. V. 2005, ApJ, 618, 874 
Herold, H. 1979, Phys. Rev. D, 19, 2868 

Herold, H, Ruder, H, & Wunner, G. 1985, Physical Review Letters, 54, 1452 
Hibschman, J. A., & Arons, J. 2001a, ApJ, 554, 624 
Hibschman, J. A., & Arons, J. 2001b, ApJ, 560, 871 
Ibrahim, A. I., et al. 2001, ApJ, 558, 237 

Johnson, M. H, & Lippmann, B. A. 1949, Physical Review , 76, 828 
Kardashev, N. S., Mitrofanov, N. S., & Novikov, I. D. 1984, AZh, 61,1113 
Kaspi, V. M., Gavriil, F. P., Chakrabarty, D., Lackey, J. R., & Muno, M. P. 

2001, ApJ, 558, 253 
Kozlenkov, A. A., & Mitrofanov, I. G. 1986, Zhurnal Eksperimental noi i 

Teoreticheskoi Fiziki, 91, 1978 
Kuiper, L., Hermsen, W., & Mendez, M. 2004, ApJ, 613, 1173 
Kuiper, L., Hermsen, W., den Hartog, P. R., & Collmar, W. 2006, ApJ, 645, 

556 

Lieb, E. H, Solovej, J. P., & Yngvason, J. 1992, Physical Review Letters, 69, 
749 



Lyutikov, M., & Thompson, C. 2005, ApJ, 634, 1223 
Lyutikov, M., & Gavriil, F. P. 2006, MNRAS, 368, 690 
Medin, Z., & Lai, D. 2007, MNRAS, 1065 

Melrose, D. B., & Parle, A. J. 1983a, Australian Journal of Physics, 36, 755 
Melrose, D. B., & Parle, A. J. 1983b, Australian Journal of Physics, 36, 799 
Mereghetti, S., Gotz, D., Mirabel, I. E, & Hurley, K. 2005, A&A, 433, L9 
Padden, W. E. P. 1994, Journal of Physics B Atomic Molecular Physics, 27, 
5419 

Potekhin, A. Y., & Pavlov, G. G. 1997, ApJ, 483, 414 
Psaltis, D., & Miller, M. C. 2002, ApJ, 578, 325 
Rea, N., et al. 2006, preprint I astro-ph/06 10802 1 

Sasaki, M., Plucinsky, P. P., Gaetz, T. J., Smith, R. K., Edgar, R. J., & Slane, 

P. O. 2004, ApJ, 617, 322 
Shabad, A. E., & Usov, V. V. 1982, Nature, 295, 215 
Shabad, A. E., & Usov, V. V. 1984, Ap&SS, 102, 327 
Shabad, A. E., & Usov, V. V. 1985, Ap&SS, 117, 309 
Shabad, A. E., & Usov, V. V. 1986, Ap&SS, 128, 377 
Silantev, N. A., & Iakovlev, D. G. 1980, Ap&SS, 71, 45 
Sokolov, A. A., & Ternov, I. M. 1968, Synchrotron Radiation (Oxford: 

Pergamon) 
Spitkovsky, A. 2006, ApJ, 648, L51 
Stix, T. H. 1976, Physical Review Letters, 36, 521 
Stumer, S. J. 1995, ApJ, 446, 292 

Sturrock, P. A. 1971, ApJ, 164, 529 

Thompson, C, ApJ, submitted I arXiv:0802.2572 1 

Thompson, C, & Blaes, O. 1998, Phys. Rev. D, 57, 3219 

Thompson, C, Duncan, R. C, Woods, P. M., Kouveliotou, C, Finger, M. H, 

& van Paradijs, J. 2000, ApJ, 543, 340 
Thompson, C, Lyutikov, M., & Kulkarni, S. R. 2002, ApJ, 574, 332 
Thompson, C, & Beloborodov, A. M. 2005, ApJ, 635, 565 
Usov, V. V. 2002, ApJ, 572, L87 
Usov, V. V., & Melrose, D. B. 1996, ApJ, 464, 306 
Wang, Z., Kaspi, V. M., & Higdon, S. J. U. 2007, ApJ, 665, 1292 
Wesson, J. 2004, Tokamaks (New York: Oxford) 
Woltjer, L. 1958, ApJ, 128, 384 

Woods, P. M., Kouveliotou, C, G6gil§, E., Finger, M. H, Swank, J., 
Markwardt, C. B., Hurley, K, & van der Klis, M. 2002, ApJ, 576, 381 

Woods, P. M., & Thompson, C. 2006, in Compact Stellar X-ray Sources, eds. 
W. H. G. Lewin and M. van der Klis (Cambridge: Cambridge University 
Press) 547 

Wunner, G., & Herold, H. 1979, Ap&SS, 63, 503 
Zhang, B., & Qiao, G. J. 1998, A&A, 338, 62 



ELECTRODYNAMICS OF M AGNETARS : PAIR CREATION AND DYNAMIC PARTICLE HEATING 2 1 





FIG. 4. — Diagrams contributing to the calculation of the cross section for 7 + 7 — > e + + e~. In case 1 (head on collision between X-ray and gamma-ray), the 
outgoing charges 1 and 2 are both in their lowest Landau state n = 0. In case 2 (gamma ray near threshold for pair single-photon pair creation), the charge that 
shares a vertex with the gamma-ray ends up in Landau state n = 0, and the charge which shares a vertex with the X-ray ends up in the state n = 1 . When the 
interaction between the two photons is resonant, there is no interference between the two channels. 



APPENDIX 

CALCULATION OF CROSS SECTION FOR 7 + 7 -» E + +E~ 

We summarize our calculation of the creation of an electron-positron pair in an ultrastrong magnetic field, via the collision of 
two photons. The basic strategy is to calculate the tree level matrix elements illustrated in Fig. |4j using magnetized wavefunctions 
for the outgoing electron and positron and the virtual electron. In this section, we use naturalized units in which fi = c= 1 . 

We focus on two special cases: 

1. - The two photons collide head on, each with a center of mass energy w\ = 0J2 > m e . Their momenta are aligned with 
the background magnetic field B. This choice of ingoing momenta corresponds approximately to the case where an energetic 
gamma-ray is created by resonant scattering of a thermal X-ray by a relativistic charge, in a magnetic field B < 4Bq. Then the 
gamma ray is below threshold for single-photon pair creation near the point of creation, u>i sin Oub.x < 2m e , but it can create a pair 
by colliding with a thermal X-ray. In the frame where the net momentum parallel to B vanishes, k\ , : + k2, z = 0, the target X-ray 
propagates nearly parallel to B, and the gamma ray propagates at a modest angle to the magnetic field. 

2. - One of the photons is close to the threshold for single-photon pair creation, 

u>\ sin^ i = 2m e - AE. (Al) 

We take this photon to propagate perpendicular to B, O^b ,1 = tt/2, and the target photon parallel to B (Fig. [TJ. The second photon 
need not have a very large energy to permit the creation of an electron and positron. 

A calculation of the cross section for 7+7 — > e + +e~ in an ultra-strong magnetic field was published by Koslenkov & Mitrofanov 
(1986) using the Johnson & Lippman (1949) electron/positron wave functions. These wave functions are not, however, eigenstates 
of the electron spin. We therefore use an alternative form derived by Sokolov & Ternov (1968) and Melrose & Parle (1983). 

The solution to the Dirac equation 

d_ 



ip = (A2) 
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is labeled by the component of the spin parallel to the magnetic field, a B = Act = ± |, the longitudinal momentum P = P ■ B, and 
the Landau level n: 

(A3) 



for the electron (positive-energy) states, and 



[^(x,t)] pna = e ip -*u%(x ± ), 



J Pmm 



for the positron (negative-energy) states. The dispersion relation is given by 

E 2 = P 2 + El=P 2 + m 2 + 2n\e\B, 



(A4) 



(A5) 



where P denotes the component of the momentum parallel to B. The wave function is localized about one of the two coordinates 
transverse to the magnetic field (the x-coordinate), 

x I = ax. 



The corresponding transverse momentum is 



P_i_ =<?xj_ x B 



-aeBy = --^y. 



In a frame where P = 0, the positive-energy wave functions are 



0*±) = 



and 







(Eq„ + m e ) l l 2 4> n {x-a) 

L(2E „ ) 1 1 1 -i(E „ -m e )^ 2 (j> n -i(x-a) 


" (Eon+nte) 1 / 2 (j) n -i{x-a) ' 
1 
L(2E „y/ 2 

. iiEon-me) 1 ' 2 (j)n(x- fl) . 

where L 3 is the normalization volume. The jc-dependence of the wave function is 



(«>0) 



(A6) 



(A7) 



(A8) 



(«>D, 



(A9) 



t>„(x-a) = 



(2"n!) 1 / 2 V Ti- 



1/4 



exp 



1 / x — a 



2 V A fi 



(A10) 



where H„ is the Hermite polynomial. In a frame boosted to a momentum P parallel to B, the wave functions become 



2L[££ () «(£ + £o„)] 1/2 



-iP{E 0n -m e ) l l 2 < 



>n-\ 



0*^ = 

The positron wave functions are given by 

v£W) = 

„C+l) 



1 

2L[EE „(E + Eo„)] 1 / 2 



1 

2L[EE „(E + E 0n )V/ 2 



(E+E 0n )(E „+m e )^ 2 (l) n 
-i(E + E „)(E „ - m e yl 2 c/>„_i 
-P(E 0n + m e ) l / 2 ^ n 

(E + E n)(E „ + m e yl 2 (j>„-\ 
-iP(E 0n -m e y/ 2 ^ 



P(E „ + m e y/ 2 ( 



i(E + E 0n )(E 0n - m e y' 1 4> n 



'n-1 
1/2, 



(All) 



(A12) 



^ n 2L[EE „(E + £„„)] V2 



i(E + E Q „)(E 0n - m e ) x l 2 c6„_i 
-P^On + m,) 1 / 2 ^ 
iP(E 0n -m e )^ 2 4> n -\ 
(E + E ()n )(E Qn + m e ) l l 2 <j) n 

P{E Qn + m e ) l l 2 <\> n - X 
-i(E+Eo„)(E 0n -m e ) 1 / 2 (f>„ 
(E + E 0n )(E 0n + m e ) 1 12 (j) n -i 

iP(E 0n -m e y/ 2 (j>„ 



(A13) 



(A14) 



In these expressions, Eo„ is the energy at P = 0, and E = E(P) (eq. HA5I ). 
The wavefunctions of the two photons (four-momenta k, = (u/;, k,), polarization vectors ef , i = 1 , 2) are normalized in the usual 

way, 
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Case 1: Head on Collision, Center-of-Momentum Frame 

We consider the photon collision not too far above the threshold for pair creation, in a range of energies where the outgoing 
electron and positron must be confined to their lowest Landau states, «_ = «+ = 0. We work in the center-of-momentum frame 
where 

u>\ = k\ z = oj; L02=—k\ z = u>. (A16) 
The energies and momenta of the outgoing electron and positron are then 



We calculate the matrix elements 



E. = E + = lu P. = -P + = ^Lo 2 -m 2 e . 

s fi =s fi m+s f m, 



(A17) 
(A18) 



S fi [l] = ie 2 J d 4 x J d 4 x' [^- 1) (x')] n =0 7 M A^(x')Gf(x'-x) 7 „Ar(x) 
S fi [2] = ie 2 J d 4 x J d 4 x' [^- 1) (x')] M =0 7^i I (^)Gf(x'-x) 7l ,A 2 -(x) 



«+=0 



«+=0 



(A19) 



using the real-space representation of the electron propagator 

G F ( x '- x )=L 2 f ^ J ^k^-o^^^xV;'!^^',*' 1 '- 1 ' 

+ i6{t - 1 ') ^ v^' ) (x')v < „';' ) (x) e , ' £ ' ( ''-' ) e -'' p(x '- x) 



(A20) 



Here ip = V>^7o = (ip*) T Jo, where * denotes complex conjugation and T the transpose. Performing the integrals over t, t' yields a 
factor 2ttS(E + + E- - uji- uj2)(E- - LJ2T EiT 1 for each of the two terms in Gf, and we have 

ie 2 1 



(A21) 
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(A22) 



These integrals involve a matrix 



where cr, are the 2 x 2 Pauli matrices. Expanding, one finds 



1 0,0 <T X \ r, Uy \ y <T t 

o -l ) £ < + u o H + u o K + u o j 



7o7^ £? = 



/ ef £ 7 \ 

£+ -e? 

£7 

\e! -a* J 



(A23) 



(A24) 



where £,- = £? ± ;'£,. 



It will be noted that only the Landau level «/ = 1 in the propagator contributes to the matrix element. That is because 

-(-1) (-1) U 7 



(A25) 
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vanishes when the two photons are incident parallel to B, and because integrals of the form 



d x±4>o(x-a{)(j) m (x-a2)e 



OC Ik; 



also vanish when m > I. We are therefore left with evaluating integrals of the form 



J 3 x</>o(x-fl 1 )^ (^-«2y (k '- p ' +P2)x = — ) 5 



2tt 



a\ — a 2 

\2 
A B 



where the result assumes k,-j_ = 0. For example, the integrals I\, I 2 take the form 
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where 
and 
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P + Pi- 
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Summing gives 



and 



r(-l)r(-l) , r(+l)r(+l)_ c .- c .+ ' 
h L 2 +1 1 L 2 - £ \ £ 2~ 



m e (P/ — P±) ( 2ir\ — / a + — a r 



A l 



2E±Ei V L 
xS(uj + P-- Pi) 5(u -P+- P,) 



\2 
A B 



2E±Ei 

x5(uj + P- + P i )S(uj-P + +P i ). 
We have used the notation P±= P- = -P+ and E± =E+ = E-, and substituted k\, z = -k 2z = u>. 
The first term in the matrix element is obtained by substituting eqs. ( IA32I ) and ( IA33t into eq. ( IA21l i. 



(2tt) 3 S(E. +E + -2u) S(P + + P.) 6 



AS 



(A26) 
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(A28) 
(A29) 
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(A31) 



(A32) 



(A33) 



(A34) 



2L 5 uj[(uj + P-) 2 +m 2 e +2eB] 

The second contribution to the matrix element is obtained by exchanging the two photons, which corresponds to taking Pj —> P--u 
in the denominator, 

„2 ,„„ _+„- 



2L 5 to[(uj-P-) 2 + m 2 e + 2eB] 



(2nY S(E. + E + - 2lu) S(P + + PJ) 6 



A, 



(A35) 



The total cross section for 7 + 7 — > e + + e , with the electron and positron both in the lowest Landau state, is obtained from 



(1 -COS#i2)<7 = 



da- 



da + f dp- f dp+ 



T J 2-K^k J 2-k\\ J 2ir J 2ir 



(A36) 



Here cos On =— 1 is the cosine of the collision angle between the two photons, and T is a normalization time. The square of a delta 
function is handled with the substitution [(2tt)5(E- + E+-2uj)] 2 -> T (2w)S(E-+ E+-2w), [{2tt)S(P++PJ}] 2 ->■ L(2ir)8(P+ + PJ), 
and so on, yielding 

pA / R \ dF.j. s+e-m 2 e\e- 2 m 2 e ' 



B 

32TTLU 2 V^Q 



dP± 



e\e 2 m 2 e 



(uj + P±) 2 +m 2 e + 2eB (oj-P±) 2 + m 2 + 2eB 

Here dE±/dP± = V± is the speed of the outgoing electron and positron. The coefficient in front is readily expressed in terms of 
the Thomson cross section 

lv( Cl ' (A38) 



(A39) 



T 3 ^ 4Tvm e 

Near threshold, the cross section vanishes if the two photons have the same linear polarization, 



a _ 3B/B Q fV± 
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Case 2: Photon 1 Near Threshold for Pair Creation 

We now turn to the case where one photon (labeled 1) is just below threshold for pair creation. We choose a frame in which 
cu\ = 2m e -AE, k\, z = 0, and L02 = ki.z- We focus on the case where there is a near resonance in the matrix element, in the sense that 
the virtual electron is nearly on mass shell. Simple kinematic considerations show that this corresponds to one outgoing particle 
in the first Landau state, with the virtual particle in level «/ = (® . The final-state particle in state n = 1 is the one that shares a 
vertex with photon 2 in Fig. [4] The final states with n_ = 1, n+ = and n_ = 0, n+ = 1 are distinguishable, and so we may calculate 
the cross section for the electron in the excited state, and then multiply by a factor 2 to obtain the total cross section. The excited 
particle also comes in two spin states a = ±1, which must be summed over. 

The matrix elements to be calculated are therefore 



Sfr^ie 2 I d 4 x I d 4 x' [^-\x')] n =llll A^x')G F (x'-x) 1 ,A , ((x) 



n+=0 



(tr_ = ±l), 



and the cross section is given by 
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where in this case cos 6*12 = 0. The condition for a near resonance is that Ej ~ E--0J2, and so we need only include the term in 
expression ( 1A20I ) for Gf(x' -x) involving the positive-energy wavefunctions. The matrix elements ( 1A40I ) become 
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where 
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Near threshold for pair creation, one has E + = Ej ~ m e , £_ ~ m e + eB/m e , and 0J2 — eB/m e . The above integrals may then be 
evaluated as 
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where 



h = e\ exp 
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and we have chosen the momentum of photon 1 to point in the y-direction. The cross section is proportional to 

E- - m e B /Bq 



| F (-1)|2 + | F (+1)|2 
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Substituting eqs. ( IA43t - (IA45b into ( IA42b . and averaging over the spins of photon 2, one finds 
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As expected, the cross section is non-vanishing only if photon 1 is polarized parallel to B. The last factor is determined by 
differentiating 

E-+E+ = (p2 + m 2 e +2eB) l/2 + (Pl + m 2 e ) 1/2 (A49) 



with P- = uj2- P+. Near threshold, 



and we have 
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The cross section can then be expressed in terms of oj (eq. l52l ). 
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